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ABSTRACT 

A  high-resolution,  multi -level  ,  primitive  equation  ocean  model  is 
used  to  examine  the  response  to  wind  forcing  of  an  idealized  flat- 
bottomed  oceanic  regime  along  an  eastern  ocean  boundary.   A  band  of 
steady  winds,  either  with  or  without  a  curl,  is  used  as  forcing  on 
both  an  f -plane  and  a  /3-plane.   In  addition,  a  stability  analysis  is 
made  to  determine  if  the  necessary  and  sufficient  conditions  for 
instability  processes  to  occur  are  satisfied.   It  is  seen  that  when 
the  wind  driven  coastal  jet  and  undercurrent  are  unstable  (which 
occurs  in  the  cases  of  wind  with  no  curl),  eddies  and  jets  are 
generated.   In  the  case  of  wind  with  curl,  since  the  Davidson  Current 
develops  rather  than  the  coastal  jet  and  undercurrent,  no  eddies 
develop.   A  comparison  of  model  results  with  available  observations 
shows  that  both  the  time -averaged  and  instantaneous  model  simulations 
of  the  coastal  jet,  undercurrent  and  eddies  are  consistent  with 
available  observational  data.   The  results  of  this  study  support  the 
hypothesis  that  wind  forcing  can  be  an  important  eddy  generation 
mechanism  for  the  California  Current  System. 
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I.   INTRODUCTION 

A.   OBJECTIVES 

The  demand  for  dependable,  accurate  ocean  prediction  models  has 
increased.  For  example,  for  U.S.  Naval  battle  group  anti-submarine 
warfare  (ASW)  tactics,  an  accurate  description  of  ocean  variability  is 
required  in  strategic  areas  of  the  world.  In  particular,  the  ability 
to  generate  spatial  and  temporal  forecasts  of  ocean  fronts,  currents 
and  eddies  is  of  crucial  importance  if  the  United  States  is  to 
maintain  superiority  in  the  ASW  arena.  One  promising  numerical  ocean 
model  has  already  been  developed  by  Miller  et  al..  (1983)  for  the  Gulf 
Stream  and  is  presently  being  used  by  the  Navy.  This  model,  which  is 
quasigeostrophic ,  can  be  initiated  with  real-time  infra-red  satellite 
and  bathythermograph  data  and  can  subsequently  develop  a  dynamical 
forecast  for  regional  Gulf  Stream  areas  (O'Brien,  1986). 

In  addition  to  the  Gulf  Stream,  the  eastern  Pacific  Ocean, 
specifically  the  California  Current  System  (CCS),  has  been  designated 
as  an  area  of  high  ASW  interest.  Models  of  the  CCS  and  other  eastern 
boundary  current  regions  are  important  not  only  for  military 
applications,  but  also  for  civilian  applications  such  as  fisheries, 
oil  recovery,  waste  disposal  and  search/rescue  operations. 

Current  numerical  modeling  efforts  of  the  CCS  at  the  Naval 
Postgraduate  School  involve  the  use  of  a  multi-level,  eddy-resolving 
primitive   equation   (PE)   ocean   model   to   investigate   generation 


mechanisms  for  the  synoptic -mesoscale  variability  of  the  CCS.  Three 
generation  mechanisms  currently  being  investigated  for  eddy  and  jet 
formations  in  the  CCS  are  wind  forcing;  topographic  forcing  due  to 
the  presence  of  the  Mendocino  Escarpment;  and  baroclinic,  barotropic 
or  mixed  instability  of  the  mean  CCS. 

The  primary  objective  of  this  thesis  is  to  investigate  the  role 
of  wind  forcing  in  the  dynamics  of  the  mean  CCS,  and  in  the  generation 
of  eddies  and  jets  in  the  region.  The  organization  of  the  thesis  is 
as  follows.  The  following  section  gives  a  brief  background  for  the 
CCS,  while  Section  II  describes  the  numerical  model  used.  The 
results  of  the  wind  forcing  experiments  along  with  a  stability 
analysis  are  presented  in  Section  III.  The  results  are  compared  with 
available  observations  in  Section  IV.  Both  a  summary  and 
recommendations  are  presented  in  Section  V. 

B.   THE  CALIFORNIA  CURRENT  SYSTEM 

The  CCS,  according  to  Hickey  (1979),  consists  of  four  currents: 
the  California  Current,  a  broad  southward  flowing  surface  current;  the 
California  Undercurrent,  a  northward  flowing  subsurface  flow;  the 
Davidson  Current,  a  surface,  poleward  flow  found  north  of  Point 
Conception;  and  the  Southern  California  Countercurrent ,  a  northward - 
flowing  surface  current  found  south  of  Point  Conception  in  the 
California  Bight.  The  ocean  circulation  along  the  California  coast  is 
forced  by  the  atmospheric  circulation  around  the  North  Pacific  High, 
the  position  and  strength  of  which  varies  seasonally  so  that  the 
offshore  pressure  gradient  along  the  west  coast  is  strongest  in  summer 
(Huyer,  1983).   These  northerly  winds  are  responsible  for  the  coastal 


upwelling  observed  from  ~April  to  October  by  driving  an  offshore 
Ekman  flux,  which  requires  a  compensating  vertical  flow  to  conserve 
mass  (Gill,  1982),  resulting  in  a  two-dimensional  (cross-shore  and 
vertical)  circulation.  Brink  (1983)  emphasized  upwelling  as  a  three- 
dimensional  process,  particularly  the  role  of  horizontal  advection  of 
heat  and  momentum  in  the  near -surface  dynamics  of  coastal  upwelling. 

Both  wind  stress  and  wind  stress  curl  are  thought  to  play 
important  roles  in  the  dynamics  of  the  mean  CCS.  At  seasonal 
frequencies,  the  alongshore  component  of  the  wind  stress  and  the  wind 
stress  curl  have  been  hypothesized  as  driving  the  seasonal  circulation 
of  the  eastern  ocean  (Carton  and  Philander,  1984).  Bryan  and  Ripa 
(1978),  Hickey  (1979),  Chelton  (1984)  and  McCreary  et  al.  (1987)  have 
stated  that  wind  stress  curl  through  Ekman  pumping  is  responsible  for 
coastal  upwelling  events  and  forces  the  pressure  field  to  establish 
associated  pressure  gradients.  Carton  and  Philander  (1984)  and 
Philander  and  Yoon  (1982)  discussed,  respectively,  how  wind  stress 
curl  and  alongshore  wind  stress  could  be  important  forcing  mechanisms 
that  contribute  to  the  seasonal  variability  in  the  CCS  near  the  coast. 
Anderson  and  Gill  (1975)  studied  the  oceanic  adjustment  resulting  from 
alongshore  wind  stress  forcing  and  concluded  that  planetary  Rossby 
waves  carry  the  coastal  upwelling  westward,  contributing  to  the  ocean 
variability.  McCreary  et  al.  (1987),  using  a  flat-bottom  ocean  model 
with  steady,  equatorward,  curl -free  wind,  showed  that  an  equatorward 
surface  coastal  jet  can  develop  due  to  the  wind  forcing  along  with  a 
poleward  undercurrent  due  to  the  pi'esence  of  an  alongshore  variation 
in  the  wind  field.   Wlien  wind  stress  curl  was  used  in  the  model  of 


McCreary  et  al_-  (1987),  a  poleward  surface  current  (i.e.,  the  Davidson 
Current)  developed  in  response  to  the  curl  near  the  coast  along  with 
an  equatorward  offshore  flow  approximately  150-175  km  from  the  coast 
(consistent  with  Chelton's  (1984)  vertical  cross -sections  of 
alongshore  geostrophic  velocity) . 

Recent  observations  have  shown  that  the  CCS  can  consist  of  intense 
meandering  current  filaments  (i.e.,  jets)  with  peak  velocities  of  ~50 
cm/sec  intermingled  with  synoptic-mesoscale  eddies  with  wavelengths  of 
several  hundred  kilometers  (Mooers  and  Robinson,  1984).  These  jets 
can  entrain  cold,  upwelled  coastal  waters  and  advect  them  hundreds  of 
kilometers  offshore.  Ikeda  and  Emery  (1984)  have  shown  that 
baroclinic  instability  associated  with  the  vertical  shear  between  the 
surface  curreiit  and  the  undercurrent  can  be  responsible  for  the  growth 
of  offshore  meander  patterns  in  the  regions  of  capes  (alongshore 
variations  in  the  coastline). 

The  role  of  wind  forcing  in  the  generation  of  eddies  and  jets  in 
the  CCS  has  not  yet  been  systematically  investigated  and  may  be  the 
most  important  generation  mechanism  for  eddy  and  jet  formations. 
Satellite  infra-red  imagery  has  shown  evidence  of  eddies  and  jets  in 
the  CCS  during  periods  of  winds  favorable  for  upwelling.  These 
observations  provide  evidence  for  wind  forcing  as  a  possible  important 
role  in  eddy  and  jet  formation.  If  wind -forced,  eddies  and  jets  could 
be  caused  by  either  a  seasonal  response  to  the  wind  field  or  short- 
lived, strong  wind  events  occurring  during  the  upwelling  season. 

Although  wind  events  may  be  important  eddy  generation  mechanisms 
(Carton,  1984;  Carton  and  Philander,  1984),  in  this  study  we  will 


examine  eddy  and  jet  formations  due  to  the  response  to  the  seasonal 
wind  field  and  leave  for  future  research  the  investigation  of  time- 
dependent  wind  forcing.  The  important  results  of  McCrearv  e_t  al  . 
(1987)  will  be  used  as  a  guide  for  the  selection  of  model  experiments. 
In  particular,  the  roles  of  both  wind  stress  and  wind  stress  curl  will 
be  examined  in  the  generation  of  eddies  and  jets  in  the  CCS. 


II.   NUMERICAL  MODEL  OF  THE  CALIFORNIA  CURRENT  SYSTEM 
A.   MODEL  DESCRIPTION 
1 .   Model  Equations 

The  numerical  model  used  in  this  research  was  developed  by 
Haney  (1985),  modified  by  Batteen  (1988a),  and  is  a  multi-level,  PE 
model.  The  model  has  ^-plane  capability,  and  uses  the  hydrostatic  and 
Boussinesq  approximations.  The  model  also  has  a  topographic 
capability,  but  in  this  study  only  a  flat-bottom  will  be  used  in  order 
to  eiisure  that  the  role  of  wind  forcing  will  be  isolated  from  the 
possible  coupled  role  of  wind  forcing  with  bottom  topography.  The 
governing  equations,  written  in  sigma  coordinates,  are  as  follows: 

a.   Momentum  Equations: 

du  -lap'  odp'dD  ^  ,       4 

dt  O^x  f^QUaoox  "' 

+     Pi;^  +      «d(")  (2-1) 


dv  _  -Idv'  odv'dD  4 

dt  "      ^o^y  '^      PoDdady      '      ^"^   '   ^"'^  ^ 


'^m' 


*  W^  +   «ci<v)  (2.2) 

b.   Continuity  Equation: 

d^d.  d\x  dv  ^  /o  T\ 

—  -t-  —  +  —  =   0  (2.3) 

do  ax  ay 


c.  Vertically   Integrated   Hydrostatic   Equation: 

po  ro         ro 

p'    =  D    I      B  de    -       I    ^[dI    Bd^]da  (2.4) 

d.  Equation  of   State: 

B  =  ag(T    -    Tq)  (2.5) 

e.  Thermodynamic   Equation: 

dT  K^a^T      Q 

—  =   -  AhV*!  +  rjT-^  +  T—    [S(z)-(w'T')] 
dt        "        D  oa      Doa 

+  5d(T)  (2.6) 

In  the  above  equations,  sigma,  denoted  by  a ,  is  equal  to  z/D. 
Equation  (2.4)  for  the  pressure  includes  the  assumption  that  the 
depth-averaged  pressure,  i.e.,  the  barotropic  mode,  is  zero.  All 
horizontal  partial  derivatives  are  on  constant  sigma  surfaces.  Table 
1  defines  the  variables  used  in  the  above  equations,  while  Table  2 
provides  other  symbols  in  the  model  equations  as  well  as  values  of 
constants  used  throughout  this  study. 
2 .   Model  Domain  and  Resolution 

The  domain  of  the  model  is  the  rectangular  region  extending 
from  124°W  to  ISO^W  and  from  36.5°N  to  42.5°N,  covering  an  area  of  6° 
longitude  by  6°  latitude  (Figure  2.1).  The  region  extends 
approximately  500  km  offshore  from  the  west  coast  of  North  America, 
and  it  spans  the  California  coastline  from  Point  Sur  in  the  south  to 
Cape  Blanco  in  the  north  (640  km).  The  horizontal  resolution  is  8  km 
in  the  east-west  direction  and  10  km  in  the  north-south  direction. 
This  horizontal  grid  resolution  should  allow  realistic,   spatial 


TABLE  1 

DEFINITIONS  OF  VARIABLES 
USED  IN  THE  MODEL 


SYMBOL  DEFINITION 

z  height  (positive  upwards) 

t  time 

p'  pressure  perturbation  from  a  vertical  average 

T  temperature 

u,v,w  eastward,  northward,  and  vertical  (sigma)  velocity 

components,  respectively 

B  buoyancy 

6^  dynamic  adjustment  term 

S  solar  radiation  in  the  ocean 

w'T'  turbulent  vertical  heat  flux 


TABLE  2 
VALUES  OF  CONSTANTS  USED  IN  THE  MODEL 


s 

% 

L 

Tq 

P 
a 

P 
o 

a 
K 

AX 
AY 
D 

<f>o 


A 


A 


0 


m 
At 

f 

g 

Am 
Ah 


M 


VALUE 


^H 
sf  c 


2  TT  day"^ 

0.958  cal  gm"^ (°K)-^ 

1.3  X  10"^ 

595.  cal  gni"^ 

278. 2°K 

1.23  X  10"^  gm  cm"^ 

1.0276  gm  cm"^ 

2.01  X  10-''(°K)'^ 

10 

8.  X  10^  cm 

1.  X  10^  cm 
4.5  X  10^  cm 
36.5°N 
42.5°N 

124. 0°W 

130. 0°W 

800.  sec 

0.93  X  lO-''  sec-^ 

980.  cm  sec'^ 

2.  X  10^^  cm''  sec-i 


2.  X  lO^''  cm"  sec"^ 


0 . 5  cm^  sec " ^ 
0 . 5  cra^  sec " ^ 
1013.25  mb 


NAME 

earth  rotation  rate 
specific  heat  of  sea  water 

drag  coefficient 
latent  heat  of  sea  water 
constant  reference  temperature 
density  of  air 
density  of  sea  water  at  Tq 
thermal  expansion  coefficient 
number  of  levels  in  vertical 
meridional  grid  spacing 
zonal  grid  spacing 
total  ocean  depth 
latitude  of  southern  boundary 
latitude  of  northern  boundary 
longitude  of  eastern  boundary 
longitude  of  western  boundary 
time  step 

Coriolis  parameter 
acceleration  of  gravity 
biharmonic  momentum  diffusion 
coefficient 

biharmonic  heat  diffusion 
coefficient 

vertical  eddy  viscosity 
vertical  eddy  conductivity 
surface  air  pressure 
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Figure    2.1  Model    domain   avea  . 
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resolution  of  mesoscale  features  in  the  CCS,  which  typically  are  on 
the  order  of  100  km  (Breaker  and  Mooers ,  1986). 

3 .  Finite  Difference  Scheme 

In  the  horizontal,  a  space -staggered  B- scheme  (Arakawa  and 
Lamb,  1977;  Batteen  and  Han,  1981)  is  used  while  a  sigma  (non- 
dimensional)  coordinate  system  is  used  in  the  vertical,  which  has  ten 
levels.  The  noise  free  version  of  the  hydrostatic  equation  in  sigma 
coordinates  advocated  by  Arakawa  and  Suarez  (1983)  and  Batteen 
(1988b)  has  been  implemented. 

4 .  Forcing  Capabilities 

The  model  can  either  be  spun  up  from  rest  by  a  surface  wind 
stress  or  heat  flux,  or  it  can  be  initialized  with  a  specified  current 
field.  In  this  study,  as  in  Renaud  (1986),  the  model  is  spun  up  using 
the  climatological  wind  fields  of  Nelson  (1977),  and  a  surface  heat 
flux  is  computed  and  used  to  damp  the  temperature  field  towards  the 
mean  climatological  value. 

5 .  Boundary  Conditions 

The  northern,  southern  and  western  boundaries  are  open  using  a 
modified  version  of  the  radiation  boundary  condition  of  Camerlengo  and 
O'Brien  (1980).  The  eastern  boundary,  representing  the  west  coast  of 
North  America,  is  closed  and  straight,  with  zero  or  free -slip 
capability . 

6 .  Heat  and  Momentum  Diffusion 

The  model  uses  biharmonic  lateral  momentum  and  heat  diffusion 
with  the  choice  of  coefficients  listed  in  Table  2.  This  formulation 
of  biharmonic,  rather  than  Laplacian  lateral  diffusion,  along  with  the 
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choice  of  coefficients  should  allow  mesoscale  eddies  to  be  generated 
via  baroclinic  and/or  barotropic  instability  processes  (Holland  and 
Batteen,  1986). 

B.   SPECIFIC  EXPERIMENTAL  CONDITIONS 
1 .   Wind  Forcing  Determination 

The  key  to  obtaining  successful  model  results  from  wind 
forcing  compared  with  CCS  observations  is  to  accurately  represent  the 
wind  field  over  the  model  domain.  Using  ship  observations  in  1 -degree 
square  areas,  Nelson  (1977)  compiled  a  complete  description  of  monthly 
wind  stress  and  wind  stress  curl  off  the  west  coast  of  the  United 
States.  These  historical  marine  wind  observations  are  utilized  in  the 
PE  model  to  determine  the  required  wind  forcing.  Since  upwelling  has 
been  observed  primarily  during  periods  of  equatorward  winds  (Mooers  e_t 
al ■ ,  1976;  Huyer ,  1983),  which  are  dominant  during  the  summer  season 
in  the  central  CCS  region,  we  will  incorporate  in  the  model  wind 
stress  and  wind  stress  curl  data  for  the  months  of  June.  July  and 
August . 

a.   Incorporation  of  Wind  Stress 

From  Figures  2.2  through  2.4,  which  show  wind  stress  data 
for  the  months  of  June  through  August,  one  can  observe  in  the  surface 
stress  field  a  large  variability  in  the  offshore  direction  and  a  small 
variability  in  the  alongshore  direction.  In  the  latter  case,  the 
winds  are  less  variable  and  parallel  to  the  coast  within  10-20%. 
Based  on  these  observations,  following  Renaud  (1986),  an  idealized 
alongshore  wind  stress  can  be  developed,  and  made  a  function  of  the 
offshore  direction.    Since  the  model  domain  has  a  straight,  north- 
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south  oriented  coastline  and  the  California  coastline  is  neither 
straight  nor  north- south  oriented  (see  Figure  2.1).  the  one -degree 
squares  of  Nelson  (1977)  were  chosen  such  that  they  were  parallel  to 
the  coast.  In  addition,  the  alongshore  component  of  Nelson's  wind 
stress  along  the  coast  was  used  as  the  model  wind  stress,  while  the 
offshore  component  was  chosen  to  be  zero  for  all  experiments. 
b.   Incorporation  of  Wind  Stress  Curl 

The  wind  stress  curl  data  for  the  months  of  Juiie  through 
August  are  shown  in  Figures  2.5  through  2.7.  A  negative  curl  is 
observed  offshore,  while  a  positive  curl  is  observed  inshore, 
resulting  in  a  line  of  zero  wind  stress  curl  parallel  to  the  coast 
approximately  200-300  km  offshore.  These  observatioiis  are  consistent 
with  Chelton  (1984),  who  observed  that  the  winds  appear  to  be 
strongest  approximately  200  km  offshore  so  that  there  is  a  nearshore 
positive  wind  stress  curl  year  round.  The  magnitude  of  this  curl 
varies  seasonally  approximatelv  in  phase  with  the  alongshore  wind 
stress . 

Following  Renaud  (1986),  wind  stress  curl  was  extracted 
from  Nelson's  climatological  data  for  the  region  of  the  FE  model 
domain.  The  wind  stress  curl  points  were  then  averaged  in  the 
alongshore  direction  and  time -averaged  for  the  three  summer  months. 
The  climatological  curve  resulting  from  the  averaging  procedures  is 
shown  in  Figure  2.8.  An  analytical,  Gaussian  function  was  then  found 
and  used  to  represent  the  averaged  curl : 

2r/ax  =  0.77exp[-0.5(XR/95.0)2]- .11,  (2.7) 
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Figure  2.2 


Climatologic^il  wind  stress  for  June  (from  Nelson 
1977).   The  inside  fiox  outlines  the  model  domain 
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Figure    2 . 3 


Climatological    wind    stress    for    Jvilv    (  t  roiii   Nelson 
1977).       The    inside    Ijox    outlines    tlif    model    cloiii.iLu 
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Figure    2.4 


Climntological    wind    stress    for    Auf;ust     (iroiii    Nelscm, 
1977).       The    inside    box    outlines    tlie    model    donuiin. 
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Figure    2 . 5 


Climatolopical    wind    stress    curl    for    June    (Iroiii    Nelson 
1977).       The    inside    l.iox    outlines    the    model    dcMiuiin. 
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Figure  2 . 6 


Cliinatological  wind  stress  curl  for  July  (from  Nelson, 
1977).   The  inside  box  outlines  the  model  domain. 
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Figure  2 . 7 


Climatologicai  wind  stress  curl  for  Avigust  (from 
Nelson.  1977).   The  inside  box  outlines  the  model 
domain . 
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where  XR  is  the  distance  in  kilometers.  The  analytical  model  curve  is 
superimposed  over  the  climatological  curve  in  Figure  2.8  and  the  two 
correlate  well . 

Since  the  PE  model  uses  wind  input  rather  than  curl  (for 
ease  of  heat  flux  and  wind  stress  calculations),  values  for  the  wind 
stress  were  obtained  at  the  temperature  gridpoints  from  integration 
offshore  of  the  analytical  curl  function  at  the  u,  v  gridpoints,  as  in 
Renaud  (1986).  Finally,  winds  were  derived  from  the  analytical  stress 
at  the  temperature  gridpoints  using  the  bulk  aerodynamic  formula. 
Plots  of  the  derived  wind  stress  and  alongshore  component  wind  (model 
versus  climatology)  are  shown  in  Figures  2.9  and  2.10,  respectively. 

c.   Incorporation  of  Alongshore  Variation  in  Wind  Stress 

Anderson  and  Gill  (1975)  and  Philander  and  Yoon  (1982) 
have  pointed  out  that  variation  in  the  alongshore  winds,  even  if  their 
curl  is  zero,  can  result  in  the  development  of  a  complex  system  of 
currents  and  contribute  to  the  seasonal  variability  of  the  CCS.  To 
examine  the  role  of  alongshore  wind  variability,  the  climatological 
alongshore  wind  stress  from  Nelson  (1977)  was  zonally  averaged  for  the 
summer  season  to  provide  the  necessary  alongshore  variation.  Figures 
2.11  and  2.12  show  plots  of  latitude  versus  wind  stress  and  alongshore 
wind  velocity,  respectively.  These  values  will  be  directly  applied  in 
two  experiments  (i.e.,  experiments  3  and  6)  as  model  stress  vice  the 
use  of  an  analytical  function  for  the  other  experiments. 
2 .   Surface  Thermal  Damping 

The  CCS  can  be  a  major  upwelling  area  during  the  spring  and 
summer  seasons.    Since  the  upwelled,  cold  water  reduces  both  the 


20 


B'O       JL'Q 


90       Q*0       ro       CO       Z'Q       VQ       00      TO- 20- 

(K>IOOT/2**I\[0/3NAa)lHn3 


Figure  2 . 8 


Wind  stress  cvirl  versus  offshore  distance-  (from 
Renaud,  U'SG). 
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Wind  stress  versus  offshore  distance  (from  Renaud, 
1986). 
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sensible  and  latent  heat  fluxes,  the  heat  flux  during  the  upwelling 
season  is  a  downward  flux  or  a  damping  factor.  As  a  result,  the  model 
ocean  temperature  fields  must  be  modified  by  surface  thermal  damping. 
Because  atmospheric  variations  in  clouds  or  air  temperature  are  not  a 
significant  forcing  (compared  to  the  wind)  mechanism  in  the  CCS, 
summer  mean  values  of  solar  insolation,  cloud  cover  and  relative 
humidity  can  be  used  in  conjunction  with  computing  sensible  and  latent 
heat  fluxes  from  bulk  formulae  (Haney  et  al..  ,  1978).  The  initial  sea 
surface  and  air  temperatures  used  were  15°C  and  IS'C,  respectively. 
Following  Renaud  (1986),  the  air  temperature  was  chosen  so  that  there 
was  not  initially  any  heat  flux  at  the  surface  when  a  steady  wind 
speed  was  used.  This  maintained  the  model  heat  budget  as  a  thermal 
damper  based  on  a  uniform  surface  wind  forcing  of  -850  cm/sec ,  and  was 
a  function  of  the  changing  sea  surface  temperature. 

Following  Haney   (1971),   a   linearized   form   of   the   model 
equations  was  used  for  the  downward  surface  flux  Q: 

Q  =  7  (T^-'^-Ts),  (2.8) 

where  T^~"  is  the  apparent  atmospheric  equilibrium  temperature,  Tg  the 
model  surf ace -layer  temperature  and  7  is  a  coupling  coefficient  which 
is  strongly  dependent  on  the  surface  wind  speed  (Haney,  1985). 
Following  Haney  (1985),  the  relaxation  time  {X'  )  for  the  surface 
temperature  can  be  related  to  the  density  of  sea  water  (p),  the 
specific  heat  of  sea  water  (Cp) ,  the  depth  of  the  mixed  layer  (h)  and 
the  winds  (7)  by  the  equation: 

A-1  ==  pCph7-^.  (2.9) 
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9   1 
Using  mean  climatological  winds  of  25  Wm'^K"-^  (Haney,   1971)  and  a 

mixed  layer  depth  of  50  m,  the  relaxation  time  value  of  around  100 

days  was  computed.   As  a  result  of  this  long  damping  time,  any  surface 

temperature  anomaly  that  may  develop  due  to  wind  forcing  should  be 

observed  long  before  complete  thermal  damping  by  the  computed  surface 

heat  flux  occurs. 

3 .   Open  Boundary  Conditions 

Numerical  models  of  coastal  ocean  regions  typically  involve 
the  treatment  of  three  open  boundary  conditions  (OBCs):  two  cross- 
shelf  and  one  offshore.  According  to  R<f>ed  and  Cooper  (1985).  the 
objective  requirements  for  a  successful,  ideal  open  boundary  condition 
are  the  following:  the  OBC  should  be  numerically  stable, 
mathematically  well-posed,  have  adequate  accuracy  and  allow 
unrestricted  fluid  motions  such  as  propagating  waves.  Furthermore, 
the  time  evolution  of  the  model's  dependent  variables  at  the  open 
boundary  and  in  the  interior  domain  should  develop  and  advect  in 
harmony,  so  that  no  disturbances  or  reflections  develop  at  the  open 
boundary . 

Chapman  (1985)  also  investigated  OBCs  using  a  barotropic, 
coastal  ocean  model  with  a  straight  coastline  and  tested  seven 
different  cross -shelf  OBCs.  He  concluded  that  there  is  no  ideal  wav 
to  handle  the  numerical  treatment  of  open  boundaries  for  models  with 
the  following  features:  two  or  three  dimensions,  rotation,  variable 
depth,  bottom  friction  or  wind  stress  forcing. 

R</>ed  and  Smedstad  (1984)  found  that  the  OBCs  of  Camerlengo  and 
O'Brien   (1980)   worked  reasonably  well   for   long,   dispersive  edge 
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waves.  For  experiments  using  alongshore  wind  stress,  the  modified 
Orlanski  radiation  OBC  of  Camerlengo  and  O'Brien  (1980)  provided  the 
best  results  in  Chapman's  (1985)  study. 

Based  on  these  previous  investigations,  we  have  chosen  to  use 
the  modified  OBCs  of  Camerlengo  and  O'Brien  (1980).  In  a  previous 
study  of  the  CCS,  Renaud  (1986)  used  the  same  OBCs  and  applied  both 
uniform  and  of f shore -varying  wind  forcing  not  only  in  the  interior, 
but  on  both  the  northern  and  southern  open  boundaries  of  the  model 
domain.  Although  an  equatorward  flowing  coastal  jet  and  poleward 
flowing  undercurrent  developed,  neither  compared  well  with 
observations  and  no  instability  of  the  mean  flow  occurred. 

In  this  study,  we  found  that,  as  in  Renaud  (1986),  if  a 
uniform  wind  forcing  is  applied  not  only  over  the  whole  domain,  but 
also  on  the  northern  and  southern  open  boundaries,  a  coastal  jet 
developed  that  was  too  deep.  We  also  found  that,  if  uniform  wind 
forcing  was  applied  everywhere  except  at  the  northern  boundary, 
coastal  Kelvin  waves  would  not  be  able  to  propagate  northward  due  to 
the  presence  of  forcing  on  the  southern  boundary.  If  uniform  wind 
forcing  was  applied  everywhere  except  at  the  southern  boundary,  the 
coastal  Kelvin  waves  would  not  be  able  to  exit  due  to  the  presence  of 
forcing  on  the  northern  boundary. 

McCreary  (1987)  showed  that  if  a  uniform  wind  stress,  as  in 
Renaud  (1986),  is  used,  a  steady  alongshore  current  will  result  that 
is  too  strong,  too  deep  and  directed  equatorward  at  all  depths.  In 
order  to  generate  a  realistic  undercurrent,  he  concluded  that  either 
an  alongshore  pressure  gradient  or  a  variation  in  the  alongshore  wind 
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stress  is  necessary.  He  recommended  and  used  wind  band  forcing,  such 
that: 

tY   =  r^YCy),  (2.10) 

where  the  meridional  profile  of  the  wind  is  confined  in  a  band  away 
from  the  northern  and  southern  boundaries.  Here  r^  is  the  wind  band 
forcing,  r^  is  the  wind  strength  for  steady  forcing,  and  Y(y)  is  the 
latitudinal  distribution  of  the  stress.  The  use  of  wind  band  forcing 
allows  an  alongshore  variation  in  the  alongshore  stress,  with  the 
result  that  a  relatively  realistic  undercurrent  can  be  set  up  bv  the 
poleward  propagation  of  coastally  trapped  waves. 

Following  McCreary  (1987),  we  have  also  imposed  a  band  of  wind 
forcing  in  the  interior  of  the  domain  in  the  north- south  direction 
away  from  both  the  northern  and  southern  boundaries.  As  we  will  see, 
when  wind  band  forcing  is  used,  both  a  realistic  coastal  jet  and 
undercurrent  can  be  generated  (compared  to  the  deep  jet  of  Renaud , 
1986). 
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III.   RESULTS  OF  WIND  FORCING  EXPERIMENTS 

A.   EXPERIMENTS  ON  AN  F- PLANE 

1 .   Experiment  I  (Uniform  Wind  Stress) 

Experiment  1  was  run  on  an  f-plane  with  a  constant  Coriolis 
parameter  of  f  q  =  9 . 3  x  10'^  sec'^  based  on  the  center  of  the  domain. 
The  wind  forcing  for  this  experiment  was  uniform,  both  in  the 
alongshore  and  cross -shore  directions  and  steady  for  a  90-day  period. 
The  constant  wind  forcing  value  of  -830.0  cm/sec  was  determined  from 
the  time-,  cross-shore-  and  alongshore -averaged  summertime, 
climatological  winds  of  Nelson  (1977),  as  described  in  Section  II. B. 
Although  the  model  had  both  free-slip  and  zero-slip  capability,  all 
experiments  had  the  eastern  boundary  condition  (west  coast  of  the 
United  States)  set  as  a  zero-slip  condition,  i.e.,  the  tangential 
velocity  was  set  equal  to  zero  at  the  coast.  This  choice  was  based  on 
the  results  of  Renaud  (1986),  who  found  that  the  zero-slip  eastern 
boundary  results  compared  more  favorably  with  observations  (Kosro, 
1987;  Huyer  and  Kosro,  1987)  than  the  free-slip  results.  The  same 
type  of  thermal  damping,  as  described  in  Section  II. B,  was  used  in  all 
experiments . 

The  initial  temperature  stratification  used  in  all  experiments 
was  a  pure  exponential  temperature  profile  as  a  function  of  height 
(z).  The  profile  (Figure  3.1)  had  a  length-scale  of  450  m  with  a 
temperature  at  the  first  sigma  level  (surface  ocean  level)  of  15°C  and 
a  value  at  great  depth  of  2°C.  The  exact  form  of  this  temperature 
function  was : 

T(z)  =  2  +  ISe^/"^".  (3.1) 
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Model  initinl  temperature  profile 
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This  temperature  profile  was  derived  from  observations  used  to 
support  the  Princeton  Dynalysis  model  (Blumberg  et.  al.,  1987)  and  is 
representative  of  the  long-term,  mean  climatological  temperature 
stratification  for  the  California  coastal  region.  In  addition,  the 
Brunt-Vaisala  frequency  profile,  N^(z),  was  calculated  analytically 
using  the  temperature  function  in  equation  (3.1)  from: 

N^  =  agdJ/dz.  (3.2) 

The  resultant  N^  profile  (Figure  3.2)  was  incorporated  for  all 
experiments . 

The  model  was  initialized  at  day  0  with  full  magnitude  winds 
and,  as  expected,  inertial  oscillations  of  near-surface  ocean  currents 
developed.  These  oscillations  soon  damped,  leaving  quasi -steady  Ekman 
offshore  transport  to  the  right  of  the  wind  stress  (Figure  3.3).  As 
explained  by  Ekman  (1905),  an  equatorward  wind  parallel  to  an  eastern 
ocean  coast  can  cause  the  surface  layer  to  move  offshore  due  to  the 
Coriolis  force.  To  maintain  continuity  of  mass,  colder  water  from 
below  the  surface  must  move  upward  to  replace  the  water  that  moves 
offshore.  As  a  result,  the  region  near  the  coast,  which  has  Ekman 
transport  as  a  major  driving  mechanism  ( Pares  -  Sierra ,  1987),  can  be 
considered  a  divergence  area. 

The  subsequent  progression  of  offshore  movement  of  cold 
temperatures,  seen  in  Figure  3.4,  is  caused  by  the  cold,  upwelled 
water,  replacing  the  offshore  movement  of  coastal  waters  due  to  Ekman 
transport.     Surface   temperature   gradients   near   the   coast   are 
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Surface  current  vectors  for  experiment  1  nt  (n) 
day  10.  (1.))  dav  20. 
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Figure  3  .  3 


(continued j  Surface  current  vectors  for  experiment  1 
iit  (c)  dav  30  and  (cl )  clay  hi). 
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~0.031°C/kin  and  slightly  increase  with  time  during  the  20-40  day 
period  shown. 

The  steady,  equatorward  wind  forcing  resulted  in  an 
equatorward,  surface  coastal  jet  (Figure  3.5)  with  a  maximum  velocity 
of  ~12-13  cm/sec  by  day  40.  This  nearshore  surface  flow  is  a  result 
of  geostrophy  and  the  thermal  wind  equations.  Gill  (1982)  showed 
that  on  an  f-plane  the  coastal  jet  should  be  confined  to  the  first 
internal  Rossby  radius  of  deformation,  which  is  ~30  km  for  the  model 
domain,  as  calculated  by  the  method  of  Feliks  (1985).  The  coastal  jet 
axis  seen  in  Figure  3.5  is  at  ~16  km,  has  a  maximum  offshore  extent  of 
~45  km  and  extends  to  ^200  m  depth.  This  coastal  jet  development 
correlates  well  with  the  steady  wind  forcing  results  of  McCrearv 
et  al.  (1987). 

A  weak,  poleward  current  of  ^1  cm/sec  is  also  seen  in  Figure 
3.5  below  the  coastal  current  at  a  depth  of  ^200- 370  m.  The  offshore 
extent  of  the  undercurrent  is  confined  to  -10  km  of  the  coast. 
McCreary  (1981)  found  that  a  poleward  undercurrent  can  develop  as  a 
result  of  either  an  alongshore  pressure  gradient  or  a  variation  in  the 
wind  stress,  and  it  can  be  established  by  the  radiation  of  Kelvin 
waves,  and  the  vertical  mixing  of  heat  and  momentum.  McCreary  (1981) 
described  the  sequence  of  events  on  the  f-plane  with  suddenly  imposed 
winds.  First,  offshore  Ekman  transport  occurred  in  the  area  of  the 
applied  wind  band.  Next,  an  upwelling  signal  propagated  rapidly 
poleward  as  a  coastal  Kelvin  wave.  As  the  Kelvin  wave  passed,  a 
coastal  jet  was  set  up  and  provided  a  source  of  water  for  the  Ekman 
transport.    Philander  and  Yoon  (1982)  found  that  Kelvin  waves  also 
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Figure    3 . 4 


Surface  contours  of  temperature  (°C)  for  experiment  1 
at  (a)  day  20.  (b)  day  30  and  (c)  day  40.   Contour 
interval  is  0  .  5°<''. . 
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Figure    3.5 


Vertical  cross-shore  sections  of  meridional  (v) 
velocity  (cm/sec)  for  experiment  1  at  (a)  day  20 
and  (b)  day  40.   Contour  interval  is  2.0  cm/sec. 
Dashed  contours  denote  southern  velocities.   The 
vertical  cross -section  was  taken  at  y  =  290  kin. 
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introduced  an  undercurrent  which  reduced  the  intensity  of  coastal 
upwelling,  but  did  not  modify  the  zonal  velocity  perpendicular  to  the 
coast . 

The  vertical  and  offshore  extent  of  the  colder,  upwelled 
waters  is  depicted  in  the  vertical  cross -sections  of  temperature 
(Figure  3.6).  The  initial  conditions  of  a  uniform,  horizontal, 
temperature  field  have  been  changed  by  the  presence  of  colder, 
upwelled  water  near  the  coast.  Consequently,  a  rise  of  isotherms, 
consistent  with  upwelling,  above  ~300  m  can  be  seen,  while  a  slight 
bending  down  of  isotherms,  consistent  with  downwelling,  below  320  m  is 
discernible.  The  near-surface  upwelled  water  extends  ~120  km 
offshore.  A  comparison  of  Figure  3.6  with  Figure  3.5  shows  that  these 
results  are  consistent  with  McCreary  (1981),  who  found  that  upwelling 
did  not  reach  deep  depths,  but  was  confined  to  above  the  core  of  the 
undercurrent;  below  the  undercurrent,  weak  downwelling  occurred. 

At  around  day  45  of  the  experiment,  the  first  evidence  of 
developing  ocean  eddies  can  be  seen  in  Figure  3.7  as  perturVjations  in 
the  zonal  current  near  the  coast  at  ~384  km.  As  will  be  seen  in 
Section  III.C,  these  eddies  develop  due  to  the  presence  of  the 
unstable,  coastal  jet.  The  time  -  sequence  of  surface  zonal  velocity, 
shown  in  Figure  3.8,  shows  the  continued  development  of  these 
mesoscale  features  in  the  center  of  the  model  domain  along  the  coast 
(between  y  ~160-448  km).  The  zonal  length  scale  maximum  for  these 
eddies  is  -'70-80  km  with  zonal  velocities  greater  than  15  cm/sec.  The 
instantaneous  zonal  velocity  field  at  day  90  (Figure  3.8c)  depicts  the 
location  of  a  weak  cyclonic  and  anticyclonic  circulation  at  y  ~320  km. 
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Figure    3 . 6 


Vertical  cross -shore  sections  of  temperature  (°C)  for 
experiment  1  at  (a)  day  20  and  (b)  day  40.   Contour 
interval  is  1.0°C.   The  vertical  cross  -  sec tion  was 
taken  at  v  =  290  l;m. 
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Figure  3.7 


Surface  contours  of  zonal  (u)  velocity  (cm/sec)  for 
experiment  1  at  day  45.   Contour  interval  i.s  2.0 
cm/sec.   Dashed  contours  denote  westward  velocities 
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Figure    3 . 8 


Surface  contours  of  zonal  (u)  velocity  (cm/sec) 
for  experiment  1  at  (a)  day  70,  (b)  day  80  and  (c) 
day  90.   Contour  interval  is  2.0  cm/sec  for  (a)  and 
(b),  5.0  cm/sec  for  (c).   Dashed  contoiu's  denote 
westward  velocities. 
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This  type  of  circulation,  known  as  a  "vortex  pair",  is  consistent  with 
CCS  observations  (Bernstein  et  al .  ,  1977). 

The  instantaneous  plots  of  meridional  velocity  plots  for  days 
70,  80  and  90  (Figure  3.9)  show  the  intensification  of  the 
equatorward,  surface  coastal  jet.  By  day  90,  the  maximum  velocity  of 
the  jet  is  greater  than  20  cm/sec  with  the  core  at  about  30  km 
offshore.   The  maximum  offshore  extent  is  ~64  km. 

The  poleward  undercurrent  (Figure  3.10)  is  shown  to  extend  to 
only  ~16-18  km  offshore.  This  undercurrent  has  a  maximum  core 
velocity  of  ~6-7  cm/sec  at  a  depth  of  ~300  m  only  ~10  km  from  the 
coast.  The  entire  undercurrent  depth  extends  from  ~60  m  to  620  m. 
The  surface  coastal  jet  is  shown  to  extend  to  ~42  km  offshore,  and 
extends  to  ~400  m  depth  offshore  of  the  undercurrent.  Offshore  of  the 
jet,  at  ~127  km,  can  be  seen  the  signature  of  an  eddy  field  at  this 
latitude.  These  eddy  fields  contribute  to  the  meandering  of  the  jet, 
as  seen  in  previous  plots. 

The  surface  dynamic  height  field,  shown  in  Figure  3.11,  was 
calculated  from  the  pressure  field  relative  to  2400  m  depth.  A 
downward  slope  of  dynamic  height  towards  the  coast  was  observed,  as 
expected.  In  addition,  an  anticyclonic  eddy  was  seen  offshore  at  50 
km  at  y  ~250  km.  A  comparison  of  Figure  3.11  with  Figure  3.12,  which 
shows  the  corresponding  temperature  fields,  shows  that  the  isotherm 
perturbations  align  favorably  with  the  generated  offshore/onshore 
geostrophic  flow  of  the  eddies.  These  results,  as  will  be  shown  in 
Section  IV,  compare  favorably  with  CCS  observations  of  cold,  intense 
seaward  flow,  called  squirts. 
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Figure    3.9 


Surface  contours  of  meridional  (v)  velocity 
(cm/sec)  for  experiment  1  at  (a)  day  70,  (I)) 
day  80  and  (c)  day  90.   Contour  interval  is  5.0 
cm/sec.   Dashed  contours  denote  southward 
velocities . 
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Figure   3 . 10 


Vertical  cross -shore  sections  of  meridional 
(v)  velocity  (cm/sec)  for  experiment  1  at  (a) 
day  70  and  (b)  day  80.   Contour  interval  is 
2.0  cm/sec.   Dashed  contours  denote  southward 
velocities.   The  vertical  cross  -  section  v^;as 
taken  at  y  =  290  km. 
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Figure    3.10 


(continued)  (c)  Vertical  cross-shore  section 
of  meridional  (v)  velocity  (cin/sec)  for 
experiment  1  at  day  90. 
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Figure    3.11 


Surface  contours  of  dynamic  height  (cm)  for 
experiment  1  at  (a)  day  70,  (b)  day  80  and  (c) 
day  90.   Contour  interval  is  2.0  cm. 
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Figure    3.12 


Surface  contours  of  temperature  (°C)  for 
experiment  1  at  (a)  day  70,  (b)  day  80  and  (c) 
day  90.   Contour  interval  is  0.5°C. 
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2.   Experiment  2  (Wind  Stress  Curl) 

Experiment  2  differed  from  experiment  1  only  in  that  the  wind 
stress  forcing  included  the  important  effect  of  a  curl.  As  discussed 
in  Section  II. B,  wind  field  data  (Figure  2.10)  derived  from  the 
analytical  representation  of  the  wind  stress  curl  (eqn.  2.7)  was  used 
to  force  the  PE  model . 

The  maximum  wind  stress  (Figure  2.10)  was  located  ~187  km 
offshore,  which,  as  expected,  corresponded  to  the  zero  value  of  wind 
stress  curl  (Figure  2.8).  Wind  stress  curl  can  cause  a  variation  in 
the  Ekman  transport  leading  to  convergence  and  divergence  of  mass, 
resulting  in  more  complex  mass  balances.  Inshore,  the  positive  wind 
stress  curl  along  the  coast  caused  Ekman  pumping  due  to  the  horizonal 
mass  divergence  at  the  surface,  which  required  water  from  below  the 
surface  to  replace  the  transported  fluid.  Offshore,  the  negative  cvu'l 
caused  horizonal  mass  convergence  which  led  to  the  occurrence  of 
downwelling . 

The  important  results  of  experiment  2  are  illustrated  in  the 
time  progression  of  surface  meridional  velocity  fields  shown  in  Figure 
3.13.  Initially,  a  surface  coastal  jet  developed  as  observed  at  day 
10  (Figure  3.13a),  but  later  (~day  40)  weakened  and  disappeared  due 
to  the  presence  of  a  dominating  poleward  current.  By  ~day  30,  a 
surface  poleward  current  was  generated  in  the  southeastern  portion  of 
the  domain.  This  current  (similar  to  the  Davidson  Current)  continued 
to  develop  during  the  entire  period,  so  that  by  day  ~90  (Figure  3.13d) 
the  current  had  reached  a  velocity  of  ~16  cm/sec  and  had  a  core 
centered  offshore  at  ~34  km.   A  broad,  slower  (~5  cm/sec),  offshore, 
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Figure    3.13 


Surface  contours  of  meridional  (v)  velocity 
(cm/sec)  for  experiment  2  at  (a)  day  10,  (b) 
day  30,  (c)  day  60  and  (d)  day  90.   Contour 
interval  is  5.0  cm/sec.   Dashed  contours 
denote  southward  velocities. 
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surface  equatorward  flow  was  discernible  at  --day  60  (Figure  3.13c)  and 
subsequently  encompassed  a  greater  extent  of  the  offshore  domain. 

The  time  sequence  of  vertical  cross-shore  sections,  shown  in 
Figure  3.14,  show  that  after  ~day  10,  the  weak  (~6  cm/sec),  surface 
coastal  jet  was  no  longer  observed.  Instead,  at  ~day  20,  a  poleward 
surface  current  developed  near  the  coast  with  its  core  located  at  a 
depth  of  ~150  m  and  a  maximum  velocity  of  ~20  cm/sec.  This  current 
extended  offshore  to  ~60  km  and  to  a  depth  of  ~1100  m.  At  ~day  60,  a 
broad,  offshore  equatorward  flow  appeared  with  a  width  of  ~174  km  and 
a  core  centered  at  ~120  km.  McCreary  et  al..  (1987)  found  similar 
results  with  his  flat-bottom  model.  In  particular,  he  found  that,  due 
to  the  presence  of  the  steady,  positive  wind  stress  curl,  both  a  deep, 
broad  poleward  surface  current  near  the  coast  and  an  equatorward 
surface  flow  located  farther  offshore  developed. 

No  eddies  or  jets  developed  in  this  experiment.  As  will  be 
discussed  in  Section  III.C,  this  is  due  to  the  lack  of  an  unstable 
coastal  surface  jet  overlying  a  poleward  undercurrent. 

3 .   Experiment  3  (Wind  with  Y- Variation  in  Wind  Stress) 

Experiment  3  was  similar  to  experiment  1  with  the  exception  of 
an  implemented  variation  in  the  alongshore  wind  stress,  as  discussed 
in  Section  II. B.  The  maximum  alongshore  wind  stress  computed  (Figure 
2.11)  was  just  south  of  the  Cape  Mendocino  region. 

The  initial  results  of  experiment  3  were  quite  similar  to 
experiment  1.  In  particular,  ocean  eddies  initially  developed  around 
days  40-50,  and  the  location  of  the  generation  (Figure  3.15)  of  these 
eddies  correlated  well  with  experiment  1.    As  seen  in  Figure  3.16, 
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Figure  3 . 14 


Vertical  cross-shore  sections  of  meridional 
(v)  velocity  (cm/sec)  for  experiment  2  at  (a) 
day  10  and  (b)  day  30.   Contour  interval  is 
2.0  cm/sec.   Dashed  contours  denote  soutliward 
velocities.   The  vertical  cross-section  was 
taken  at  y  =  290  km. 
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(continued)  Vertical  cross  -  shore  sections  of 
meridional  (v)  velocity  (cm/sec)  for 
experiment  2  at  (c)  day  60  and  (d)  day  90. 
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which  shows  surface  contours  of  zonal  and  meridional  velocity,  dynamic 
height  and  temperature  fields  at  day  90,  the  core  of  the  equatorward 
coastal  jet  is  centered  ~32  km  offshore  with  a  maximum  velocity  of  ~20 
cm/sec . 

A  notable  difference  from  experiment  1  was  that  experiment  3 
did  not  have  the  eddy  development  occur  as  far  north  as  in  experiment 
1,  i.e.,  no  eddies  developed  north  of  ~385  km.  This  was  likely  due 
to  the  maximum  wind  stress  occurring  south  of  this  region. 
Temperature  perturbations  associated  with  offshore  jets  and  eddies 
also,  as  expected,  occurred  further  to  the  south.  By  ~day  45,  the 
surface  coastal  jet  (Figure  3.17)  had  intensified  to  ~12  cm/sec  with 
the  core  ~25  km  offshore.  The  maximum  offshore  extent  of  the  jet  was 
~63  km.  A  weak  poleward  undercurrent  (~2  cm/sec)  was  seen  at  a  depth 
of  -200-400  m.  The  meridional  velocity  field  (Figure  3.17b)  at  day  90 
showed  that  the  undercurrent  velocity  had  increased  slightly  to  ~4 
cm/sec  while  the  surface  coastal  jet  velocity  had  increased  to  ~14 
cm/sec.  A  weak  surface  poleward  velocity  was  also  observed  at  ~87  km 
offshore  of  the  jet.  which  was  associated  with  the  development  of  an 
ant icy clonic  eddy. 

This  experiment  showed  that  the  variation  in  alongshore  wind 
stress  can  play  a  key  role  in  the  location  of  the  eddy  development 
region.  In  particular,  eddies  can  be  generated  in  the  localized  area 
of  maximum  alongshore  wind  stress. 
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Figure  3.15 


Surface  contours  of  zonal  (u)  velocity 
(cm/sec)  for  experiment  3  at  day  50.   Contour 
interval  is  2.0  cm/sec.   Dashed  contoiu's 
denote  westward  velocities. 
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Figure    3.16 


Surface  contours  of  (a)  zonal  (u)  velocity 
(cm/sec),  (b)  meridional  (v)  velocity 
(cm/sec),  (c)  dynamic  height  (cm)  and  (d) 
temperature  (°C)  for  experiment  3  at  day  90. 
Contour  interval  for  (a)  aiid  (b)  is  5.0 
cm/sec.   Contour  interval  for  (c)  is  2.0  cm 
and  (d)  is  0.5°C.   Dashed  contours  denote 
westward  velocities  in  (a),  southward 
velocities  in  (b)  and  negative  values  in  (c) 
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Figure    3.17 


Vertical  cross-shore  sections  of  meridional 
(v)  velocity  (cm/sec)  for  experiment  3  at  (a) 
day  45  and  (b)  day  90.   Contour  intei'val  is 
2.0  cm/sec.   Dashed  contours  denote  southward 
velocities.   The  vertical  cross  -  sec tion  was 
taken  at  y  =  290  kin. 
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B.   EXPERIMENTS  ON  A  BETA- PLANE 

1 .   Experiment  4  (Uniform  Wind  Stress) 

Experiment  4  used  the  same  parameters  and  forcing  mechanisms 
as  in  experiment  1,  but  used  a  ^-plane  rather  than  an  f -plane.  The 
Beta  effect,  i.e.,  the  variation  of  the  Coriolis  parameter  with 
latitude  due  to  the  curvature  of  the  earth,  allows  the  existence  of 
freely  propagating  planetary  waves,  i.e.,  Rossby  waves  (Gill,  1982). 
Due  to  the  Beta  effect,  the  surface  coastal  jet  is  not  required  to  be 
confined  within  a  Rossby  radius  of  deformation  of  the  coast  (McCreary 
et  al  ■  .  1987).  The  offshore  radiation  of  Rossby  waves,  according  to 
McCreary  (1987),  can  contribute  to  the  generation  of  an  alongshore 
pressure  gradient  field,  which,  as  mentioned  previously,  can  cause  a 
poleward  undercvirrent  to  develop. 

The  results  of  experiment  4  showed  the  generation  of  eddies  at 
~day  40,  as  seen  in  the  zonal  velocity  plots  of  Figure  3.18.  The  eddy 
perturbations  were  generated  a  little  farther  to  the  north  on  the  /9- 
plane  than  in  experiment  1.  Zonal  eddy  velocities  ranged  from  ~5-12 
cm/sec,  at  ~day  70  (Figure  3.18d).  The  time  progression  of  the 
westward  propagation  of  the  faster,  longer  wavelength  Rossby  waves  can 
be  seen  in  Figures  3.18a-d.  As  expected,  progressively  shorter 
wavelength  Rossby  waves  propagate  offshore  in  subsequent  times. 
Surface  contour  plots  for  instantaneous  velocity,  temperature  and 
dynamic  height  fields  at  day  90  are  shown  in  Figure  3.19.  The  maximum 
zonal  velocity  (Figure  3.19a)  reached  is  ~15  cm/sec  coinciding  with 
the  intense  generation  of  an  anticyclonic  eddy  at  y~224  km  (Figure 
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Figure    3.18 


Surface  contours  of  zonal  (u)  velocity 
(cm/sec)  for  experiment  4  at  (a)  day  40,  (h) 
day  50,  (c)  day  60  and  (d)  day  70.   Contour 
interval  is  2.0  cm/sec.   Dashed  contours 
denote  westward  velocities. 
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3.19c).    The  13.5°C  temperature  anomaly  in  the  same  area   (Figure 
3.19d)  is  likely  associated  with  the  warm  core  of  this  anticyclone. 

The  surface  coastal  jet  with  equatorward  flow  (Figure  3.19b) 
has  maximum  velocities  of  ~10-15  cm/sec  which  are  located  offshore  at 
~56  km.  This  was  approximately  25  km  farther  offshore  and  5  cm/sec 
weaker  than  the  results  of  experiment  1.  A  comparison  of  the  coastal 
jet  in  this  experiment  at  day  90  (Figure  3.20)  and  that  in  Experiment 
1  (Figure  3.10c)  shows  that  the  coastal  jet  is  shallower  and  the 
undercurrent  is  deeper  in  this  experiment.  The  maximum  offshore 
extent  of  the  undercurrent  is  ~30  km  ,  which  is  -10  km  farther 
offshore  than  in  experiment  1. 

This  experiment  showed  that,  due  to  the  Beta  effect, 
circulation  does  not  necessarily  have  to  be  confined  within  a  Rossby 
radius  of  deformation  of  the  coast.  As  in  previous  no  curl 
experiments,  a  coastal  jet  develops  within  a  Rossby  radius  of 
deformation  of  the  coast  (which  is  ~30  km,  based  on  the  method  of 
Feliks,  1985),  within  which  eddies  develop.  Due  to  the  use  of 
vertical  mixing  in  the  model,  the  Rossby  waves  do  not  advect  the  jet 
or  eddies  offshore.  Rather,  the  Rossby  waves  propagate  offshore  as  a 
packet,  with  longer  and  then  shorter  waves  moving  westward,  as 
expected . 

2 .   Experiment  5  (Wind  Stress  Curl) 

This  experiment  paralleled  experiment  2,  but  used  a  /3-plane 
rather  than  an  f -plane.  Surface  contours  of  zonal  velocity,  shown  in 
Figure  3.21,  depict  the  offshore  radiation  of  Rossby  waves,  as 
expected  with  the  inclusion  of  the  beta  effect,   and  the  initial 
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Figure    3.19 


Surface  contours  of  (a)  zonal  (u)  velocity 
(cm/sec),  (b)  meridional  (v)  velocity 
(cm/sec),  (c)  dynamic  height  (cm)  and  (d) 
temperature  (°C)  for  experiment  4  at  day  90. 
Contour  interval  for  (a)  and  (b)  is  5.0 
cm/sec.   Contour  interval  for  (c)  is  2.0  cm 
and  (d)  is  0.5°C.   Dashed  contours  denote 
westward  velocities  in  (a),  southward 
velocities  in  (b)  and  negative  values  in  (c) 
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Figure  3.20 


Vertical  cross -shore  sections  of  meridional 
(v)  velocity  (cm/sec)  for  experiment  4  at  day 
90.   Contour  interval  is  2.0  cm/sec.   Dashed 
contours  denote  southward  velocities.   The 
vertical  cross-section  was  taken  at  y  =  290 
km. 
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generation  of  an  eddy  field  between  y~500-575  km.  This  experiment 
differed  from  experiment  2  (wind  stress  curl  on  an  f -plane)  which  had 
neither  of  these  features.  It  also  differed  from  experiment  1.  (no 
curl  on  an  f- plane)  which  had  of  course  no  Rossby  waves,  but  had 
eddies  developing  much  earlier  ('-day  45  rather  than  '-day  80)  and 
farther  south  (~100  km).  The  generation  of  the  eddies  here,  as  in 
previous  experiments  in  which  eddies  were  generated,  was  due  to  the 
presence  of  the  unstable  coastal  jet  and  undercurrent,  as  will  be 
shown  in  Section  III.  C. 

A  time  progression  of  surface  contours  of  meridional 
velocity,  shown  in  Figure  3.22,  showed  the  weak  surface  coastal  jet 
confined  to  within  '-16  km  from  the  coast  in  the  northern  portion  of 
the  domain.  A  poleward  current  surfaced  at  '-day  60  and  extended  from 
the  northern  to  the  southern  domain  boundaries  by  ~day  90,  with  a 
maximum  velocity  of  ~10-12  cm/sec.  This  current  was  weaker  ('-5 
cm/sec)  and  extended  ~16  km  farther  offshore  than  in  experiment  2(f- 
plane ) . 

The  broad  equatorward  flow  at  '-y  =  100-200  km  offshore  was 
best  observed  in  the  vertical  cross -shore  sections  of  meridional 
velocity  (Figure  3.23).  The  flow  developed  by  -'day  30  and  reached  a 
maximum  velocity  of  ~5  cm/sec  by  '-day/90.  The  stronger  surface 
coastal  jet  weakened  throughout  the  time  period  and  became  a  weaker, 
near-coastal  feature  as  the  poleward  undercurrent  grew  and  surfaced. 
The  core  of  the  poleward  current  was  '-8  cm/sec  weaker  and  '-15  km  wider 
than  in  experiment  2 . 
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Figure    3.21 


Surface  contours  of  zonal  (u)  velocity 
(cm/sec)  for  experiment  5  at  (a)  day  40  and 
(b)  day  80.   Contour  interval  is  2.0  cm/sec 
Dashed  contours  denote  westward  velocities. 
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Figure    3.22 


Surface  contours  of  meridional  (v)  velocity 
(cm/sec)  for  experiment  5  at  (a)  day  10,  (1)) 
day  30,  (c)  day  60  and  (d)  day  90.   Contour 
interval  is  5.0  cm/sec.   Dashed  contours 
denote  southward  velocities. 
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Figure   3.23 


Vertical  cross -shore  sections  of  meridional 
(v)  velocity  (cm/sec)  for  experiment  5  at  (a) 
day  10  and  (b)  day  30.   Contour  interval  is 
2.0  cm/sec.   Dashed  contours  denote  southward 
velocities.   The  vertical  cross  -  section  is 
taken  at  y  =  550  km. 
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Figure    3.23 


(continued)  Vertical  cross-shore  sections  of 
meridional  (v)  velocity  (cm/sec)  for 
experiment  5  at(c)  day  60  and  (d)  day  \K) . 
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The  offshore  side  of  this  current  also  had  a  much  weaker 
vertical  shear  than  that  of  experiment  2;  however,  unlike  experiment 

2,  due  to  the  presence  of  the  coastal  jet,  eddies  were  generated. 

This  experiment  showed  that,  due  to  the  presence  of  the  curl, 
both  the  Davidson  Current  and  an  equatorward  flow  more  than  100  km 
offshore  can  develop,  as  in  experiment  2.  Due  to  the  inclusion  of  the 
/3-plane  rather  than  the  f -plane  (as  in  experiment  2),  a  coastal  jet 
inshore  of  the  Davidson  Current  also  develops,  and  within  this  jet, 
eddies  are  generated.  Due  to  the  Beta  effect,  eddies  develop  further 
north  than  in  experiment  1,  and,  as  expected,  Rossby  waves  propagate 
offshore  as  a  packet,  with  longer  waves  followed  by  shorter  waves. 
This  experiment,  because  of  the  inclusion  of  both  wind  stress  curl  and 
the  /9-plane,  should  be  more  representative  of  observations  in  the  CCS, 
than  in  the  other  experiments . 

3 .   Experiment  6  (Wind  with  Y-Variation  in  Wind  Stress) 

This  experiment  used  the  same  forcing  parameters  as  experiment 

3,  but  used  a  ^-plane  rather  than  an  f -plane.  In  general,  the  results 
(Figures  3.24,  3.25)  were  similar  to  experiment  3;  however,  with  the 
/9-plane,  eddies,  as  seen  in  the  zonal  velocity  fields  at  day  50 
(Figure  3.24)  developed  ~20-25  km  farther  north,  as  expected.  The 
surface  coastal  jet,  seen  in  Figure  3.25,  had  a  maximum  core  velocity 
of  ^15  cm/sec  and  was  centered  offshore  ~32  km  in  the  north  and  -'64  km 
in  the  central  part  of  the  model  domain.  The  maximum  offshore  extent 
of  the  surface  jet  was  ~16  km  farther  offshore  than  in  experiment  3. 
The  vertical  cross -shore  section  of  meridional  velocity  (Figure  3.26) 
showed  the  time  progression  of  both  the  coastal  jet  and  the  poleward 
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Figure    3.24 


Surface  contours  of  zonal  (u)  velocity 
(cm/sec)  for  experiment  6  at  day  50.   Contour 
interval  is  2.0  cm/sec.   Dashed  contours 
denote  westward  velocities. 
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Figure    3.25 


Surface  contours  of  (a)  zonal  (u)  velocity 
(cm/sec),  (b)  meridional  (v)  velocity  (cm/sec), 
(c)  dynamic  height  (cm)  and  (d)  temperature  (°C) 
for  experiment  6  at  day  90.   Contour  interval  for 
(a)  and  (b)  is  5.0  cm/sec.   Contour  interval  for 
(c)  is  2.0  cm  and  (d)  is  0.5°C.   Dashed  contours 
denote  westward  velocities  in  (a),  southward 
velocities  in  (b)  and  negative  values  in  (c). 
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Figure    3.26 


Vertical  cross-shore  sections  of  meridional  (v) 
velocity  (cm/sec)  for  experiment  6  at  (a)  day  45 
and  (b)  day  90.   Contour  interval  is  2.0  cm/sec. 
Dashed  contours  denote  southward  velocities.   The 
vertical  cross -section  was  taken  at  y  =  290. 
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undercurrent.  Both  currents  were  weaker  than  their  counterparts  in 
experiments  1  and  3.  The  core  of  the  undercurrent  was  also  deeper  and 
farther  offshore  than  in  the  other  experiments,  consistent  with  the  (i 
effect . 

This  experiment  showed  that  both  the  Beta  effect  and  the 
variation  in  alongshore  wind  stress  can  play  key  roles  in  the  location 
of  the  eddy  development  region.  When  the  variation  in  alongshore  wind 
stress  is  included,  eddies  are  generated  in  the  localized  area  of 
maximum  alongshore  wind  stress,  which  here  is  south  of  Cape  Mendocino. 
With  the  addition  of  the  Beta  effect,  eddies  are  also  generated  north 
of  this  area. 

C.   STABILITY  ANALYSIS 

Emery  and  Mysak  (1980)  and  Thompson  (1984)  provided  support  for 
baroclinic  instability  as  an  important  generation  mechanism  for  eddies 
off  Vancouver  Island.  Ikeda  and  Emery  (1984)  used  satellite  imagery 
off  Oregon  and  Northern  California  to  show  the  growth  of  coastline, 
cape  -  initiated ,  current  meanders  as  a  result  of  baroclinic  instability 
of  the  surface  coastal  jet  and  California  undercurrent.  These 
meanders  developed,  moved  offshore,  were  cutoff  and  eventually  became 
pairs  of  isolated  eddies  or  "vortex  pairs"  (Bernstein  §X  a\.,  1977). 
Mooers  and  Robinson  (1984)  showed  that  an  intense  offshore 
quasigeostrophic  jet  was  often  located  between  these  cyclonic  eddy- 
anticyclonic  eddy  pairs  that  entrained  cold,  upwelled  water  and 
rapidly  (50  cm/sec)  advected  it  hundreds  of  kilometers  offshore. 

Evidence  for  both  barotropic  and  baroclinic  processes  as  important 
generation  mechanisms  for  eddies  was  provided  by  Thompson  (1984),  who 
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showed  that  the  contribution  of  both  instabilities  was  required  for 
the  generation  of  a  cyclonic  eddy  off  Vancouver  Island.  Batteen 
(1988a)  used  a  simulated  CCS  surface  coastal  jet  and  undercurrent  and 
showed  that  the  generation  of  complex  eddy  and  jet  patterns  could  be 
attributed  to  the  baroclinic  instability  created  by  the  vertical  shear 
between  the  two  currents. 

Despite  these  Investigations,  the  dynamical  reasons  for  the 
generation  of  the  complex  eddy  and  jet  patterns  of  the  CCS  is  not 
completely  understood,  and  so  will  be  examined  here.  First,  we 
examine  the  necessary  conditions  to  determine  the  potential  for  the 
flow  field  to  become  unstable  and  develop  eddies  and  jets.  It  is 
known  that  a  necessary  condition  for  barotropic  instability  on  an  f- 
plane  is  a  sign  change  in  the  velocity  profile  curvature  (Haltiner  and 
Williams,  1980).  Since,  in  all  the  experiments  of  this  study,  a 
maximum  velocity  core  was  located  offshore,  either  due  to  the 
equatorward,  surface  coastal  jet  or  the  poleward  undercurrent  (that 
surfaced  in  experiments  2  and  5).  the  necessary  condition  for 
barotropic  instability  to  occur  was  satisfied. 

In  addition  to  barotropic  instability,  baroclinic  instability 
could  also  be  significant  due  to  the  available  energy  from  the 
vertical  shear  of  the  coastal  jet  and  undercurrent.  Watts  (1983) 
examined  the  distributions  of  potential  vorticity  in  the  Gulf  Stream 
as  a  signature  of  baroclinic  instability.  The  cross-stream 
distribution  of  potential  vorticity,  which  is  a  conservative  quantity 
(Watts,  1983),  was  used  to  evaluate  the  necessary  conditions  for 
instability  associated  with  quasigeostrophic  theory.    Gill   (1982) 
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stated  that  without  these  necessary  conditions,  the  required  available 
energy  would  not  be  released  so  that  instability  would  not  develop.  A 
necessary  condition  for  baroclinic  instability  to  occur  is  that  the 
cross- stream  derivative  of  potential  vorticity  must  change  sign 
somewhere  within  the  domain.  In  addition,  the  product  of  the  cross- 
stream  derivative  and  the  basic  current  are  required  to  be  positive. 
Finally,  the  coastal  jet  must  meet  the  requirement  for  linearization 
of  a  basic  state  current  that  is  slowly  changing  in  space  and  time 
(Robinson,  1983);  this  requirement  is  met  by  the  structure  of  the 
coastal  jet  in  this  study. 

Watts  (1983)  examined  the  potential  vorticity  (q)  signature  in 
the  Gulf  Stream  using  the  following  expression,  in  cartesian 
coordinates : 

/ir     ,.  \  ^     3T  3V  /TON 


where 


av     au  /o  / N 

C  =  c,   -  a  (3-4) 

ax        ay 


Following  Watts,  we  similarly  studied  the  coastal  jet  to  determine  its 
potential  for  instability.  A  cross- section  plot  of  time -averaged 
(days  30-40)  and  meridionally  averaged  potential  vorticity  (Figure 
3.27)  for  experiment  1  showed  the  tendency  for  potential  vorticity  to 
be  uniform  along  isothermal  surfaces  and  to  change  vertically, 
consistent  with  the  offshore  temperature  stratification.  The  time- 
average  of  days  30-40  was  chosen  because  it  was  the  period  during 
which  the  instability  occurred.  The  range  of  the  potential  vorticity 
was  between  0.0-2.2  x  10^°C  m'^  s"^  offshore  of  the  coastal  jet.  A 
relative  minimum  existed  in  the  surface  layer  of  the  offshore  region 
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Figure  3.27     Vertical  cross  -  section  of  potential  vorticity 

(°C  m'^  s'M  for  the  time-averaged  days  30-40  of 
experiment  1,  scaled  by  lO'' .   Contour  interval  is 


0.2  °C  m 
taken  at 
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The  vertical  cross  -  section  was 


290  km. 
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due  to  weak  stratification  from  turbulent  vertical  mixing.  Strong 
upwelling  in  the  nearshore  region  caused  weak  stratification  and 
deeper  minimum  values.  A  relative  maximum  of  potential  vorticity  was 
located  at  a  depth  of  ~100  m,  greater  than  63  km  offshore,  which 
corresponded  to  the  "seasonal"  thermocline  in  the  model.  (All 
experiments  in  this  study,  where  instability  occurred,  showed  similar 
potential  vorticity  results,  and  so  will  not  be  shown  here). 

The  cross-stream  derivative  of  potential  vorticity  was  plotted 
(Figure  3.28)  to  examine  the  necessary  condition  for  baroclinic 
instability  that  the  cross-stream  derivative  of  potential  vorticity 
must  change  sign  somewhere  in  the  domain.  The  plot  was  obtained  by 
first  calculating  the  horizontal  gradient  (9q/5x)  and  then  multiplying 
by  one  grid  length  (Ax)  for  scaling  purposes.  From  Figure  3.28,  it  is 
obvious  that  the  cross-stream  derivative  meets  the  necessary  condition 
(of  a  sign  change)  for  instability  in  the  vicinity  of  the  coastal  jet 
and  undercurrent. 

Renaud ' s  (1986)  potential  vorticity  analysis  had  similar  results 
as  this  study;  however,  even  with  both  barotropic  and  baroclinic 
necessary  conditions  met,  his  experiments  did  not  show  any  development 
of  instability.  The  boundary  conditions  of  our  study,  unlike  Renaud 
(1986),  used  a  band  of  wind  forcing  in  the  middle  of  the  domain,  as 
described  in  Section  II. B. 3,  which  was  important  in  the  generation  of 
a  more  realistic  coastal  jet  and  undercurrent,  which  resulted  in  the 
generation  of  mesoscale  eddies  and  jets. 

Following  Renaud  (1986),  a  simple  baroclinic  instability  model 
(Holton,   1979)  was  used  to  examine  the  sufficient  conditions  for 


76 


3     S 


CM 
CN 


CM 

to 


CM 

to 

5 


CM 

to 


(uj)  ^]d^Q 

CM     CM     CM 
to     to     to 

5   F^    CO 


CM 

to 

CM 

CM 

CM 

CM 

CM 

to 
o 

to 

to 

to 

5 

to 

M; 

tM 

a) 

o 

c: 

o 

o-j 

to 

X) 

CM 

-t-f 

V) 

u 

1 

1 
-t-J 

in 

^ 

CM 

> 

f) 

Figure  3.28    Vertical  cross-section  of  the  cross -stream 

derivative  of  potential  vorticity  (°C  m"^  s"')  for 
the  time -averaged  days  30-40  of  experiment  1, 
scaled  by  10^.   Contour  interval  is  0.1  °C   m"'  s'' 
Dashed  contours  denote  negative  values.   The 
vertical  cross  -  sec  tion  was  tal:en  at  v  =-  2^'0  km. 
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baroclinic  instability.  The  ocean  was  divided  into  two  discrete 
layers.  The  upper  layer  extended  from  the  surface  to  150  m ,  which 
coincided  with  the  bottom  of  the  coastal  jet,  while  the  lower  layer 
extended  from  150  m  to  300  m,  basically  to  the  maximum  velocity  core 
of  the  poleward  undercurrent.  The  arrangement  of  variables  in  the 
vertical  for  the  two -layer  baroclinic  model  is  shown  in  Table  3. 

TABLE  3 

ARRANGEMENT  OF  VARIABLES  IN  THE  VERTICAL 
FOR  THE  TWO -LAYER  BAROCLINIC  MODEL 

0  m   wq=0  level  0 

-75  m   Vl  level  1 

-150  m  W2  level  2 

-225  m  ._ V'S  level  3 

-300  m  W4=0  level  4 

The  quasigeostrophic  (QG)  vorticity  equation  was  applied  at  leA'^els 
1  and  3,  whereas  the  QG  thermodynamic  equation  was  applied  at  level  2. 
The  streamf  unction ,  V'2  "^^  computed  by  linear  interpolation  between 
levels  1  and  3.  This  completed  a  closed  set  of  prediction  equations 
with  /9=0 ,  as  in  experiments  1  through  3,  using  the  PE  model  of  this 
s  tudy . 

The  length  scale  used  in  this  model  was  the  Rossby  radius  of 
deformation  (\'^)    defined  by: 

f  ^ 

where  Az  was  the  layer  thickness  (150  m)  and  N^  was  the  stability 
parameter  defined  by  equation  (3-2). 
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Following  Renaud  (1986),  the  phase  speed  (c)  was  computed  by 
solving  the  closed  set  of  QG  equations,  by  assuming  wave  -  type 
solutions,  thus 

c  =  Vm±  Vt  [^^]    ^/^  (3-6) 

where  the  vertically  averaged  current  (Vj^)  was  defined  as: 

V   ,  V 
1+3 
V^  = (3-7) 

and  the  basic  state  "thermal  wind"  (Vj)  was  defined  as: 

^  +  ^3 
Vt  =  ^ (3.8) 

For  waves  with  meridional  wavenumbers  (k)  that  satisfied  k^  <  2A^. 

there  was  an  imaginary  component  of  the  phase  speed  (C^).   Based  on 

this  criterion,  the  critical  wavelength  (L^)  was  defined  as  follows: 

Lc  =  2  ^/"(tt/A)  ,  (3-9) 

or 

L^  =  [(Az)7r(2N2)^/2]/fo  .  (3-10) 

For  the  experiments  on  an  f -plane  (^=0),  the  criteria  for  instability 
did  not  depend  on  the  magnitude  of  the  basic  state  "thermal  wind" 
(Vj) ;  therefore  all  wavelengths  longer  than  Lq  with  even  the  slightest 
vertical  shear  were  unstable. 

Any  wave  that  was  longer  than  this  critical  wavelength  amplified 
at  the  growth  rate  (6)   where: 

5  =  k  Ci  .  (3-11) 
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The  growth  rate  (6)  depended  on  V^  and  N^,  was  a  function  of  the 
meridional  wavenumber  (k)  and  was  a  maximum  at  an  intermediate 
wavenumber  (0<k'^<2A^).  The  e-folding  growth  time  was  computed  by  the 
inverse  of  8 . 

Batteen  et.  al.  (1988)  investigated  the  generation  of  eddies  by 
baroclinic  instability  due  to  an  imposed  initial  baroclinic  jet,  that 
represented  the  mean  CCS  during  the  upwelling  season.  A  comparison  of 
this  unstable  jet  (BCJET)  was  made  with  experiment  1  to  determine  the 
similarities  or  differences  with  a  known  baroclinic  instability 
example.  For  the  established  unstable  jet,  computations  were  made 
using  the  time  -  (days  1-10),  zonal  -  (within  64  km  of  the  coast)  and 
meridional  -  averaged  data  for  two  layers  (0-250  m  and  250-500  m)  . 
The  same  procedure  was  used  for  experiment  one;  however  we  time- 
averaged  for  days  1-10  (the  initial  establishment  of  the  coastal  jet 
overlying  the  poleward  undercurrent)  and  days  30-40  (the  period 
immediately  preceding  the  generation  of  eddies).  Additionally,  since 
the  mean  currents  were  shallower  and  did  not  extend  as  far  offshore, 
we  used  the  two  layers  described  earlier  of  0-150  m  and  150-300  m  and 
varied  the  zonally  averaged  data  to  correspond  to  the  offshore  extent 
of  the  mean  currents  (~8-16  km).  The  ten-day  averaged  data  was  then 
input  to  the  two  layer  baroclinic  model  to  determine  the  baroclinicity 
of  the  flat  bottom,  f-plane  experiments  with  a  coastal  jet  and 
undercurrent . 

Table  4  shows  the  internal  Rossby  radius  of  deformation  {X'^),  the 
critical  wavelength  (Lq)  ,  the  basic  state  "thermal  wind"  (V-j-)  and  the 
most  unstable  wavelength  (h^)   with  its  associated  e-folding  time  (6'^) 
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for  the  unstable  jet  and  experiment  1.  A  comparison  of  the  wind 
forcing  results  with  that  of  the  unstable  jet  in  Table  1  show  that  the 
results  are  comparable.  The  short  e-folding  times  show  that  once  this 
unstable  state  has  been  reached,  eddy  development  should  occur  after 
about  a  week  of  model  integration. 

The  results  of  experiment  1  were  as  expected  from  Table  1;  that 
is,  as  the  vertical  shear  increased  with  time  due  to  the  developing 
coastal  jet  and  undercurrent,  the  shorter  waves  became  the  more 
unstable  waves  with  shorter  e-folding  times.  However,  it  should  be 
noted  that  these  e-folding  times  are  only  an  approximation  to  the 
actual  baroclinic  growth  rates  since  a  simple  two-layer  baroclinic 
model  \<!as    used  for  the  computations. 

TABLE  4 
TWO -LAYER  BAROCLINIC  MODEL  RESULTS 
exp.(days)      A'^km)      L^^Ckm)     Vj(cm/s)     Ljj^(km)     6"' (days) 


BCJET(l-lO) 

9.60 

42.7 

9.3 

67.7 

2.05 

1(1-10) 

7.12 

31.6 

3.0 

51.6 

4.70 

1(30-40) 

6.50 

28.9 

3.7 

43.9 

3.45 

The  baroclinic  instability  growth  rate  of  the  unstable  jet  is 
shown  in  Figure  3.29  as  a  function  of  meridional  wavelength  (L=27r/K)  . 
The  shorter  wavelengths  had  the  fastest  growth  rate  compared  to  the 
longer  wavelengths;  however  the  difference  in  their  growth  rates  was 
not  that  significant.   The  preferred  growth  rate  occurred  at  the  most 
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unstable  wavelength  of  ~70  km,  which  correlates  well  with  Table  4. 
The  plot  of  e-folding  time  versus  meridional  wavelength  (Figure  3.30) 
for  experiment  1  of  the  wind  forcing  results  showed  the  time 
progression  of  the  e-folding  rates  with  the  development  of  the  coastal 
jet  and  undercurrent.  The  large  difference  between  days  1-10  and  days 
11-20  was  due  to  the  time  for  the  wind  field  to  initially  set  up  the 
mean  current  features  (instead  of  being  imposed  in  the  model  as  in 
Batteen  et,  al..  ,  1988).  At  the  shorter  wavelengths  (below  ~80  km)  the 
e-folding  times  are  similar  for  the  days  11-50,  with  the  most  unstable 
wavelength  of  ~45  km  corresponding  to  a  growth  time  of  ~3 . 5  days.  The 
longer  wavelength  shows  a  small  growth  time  difference  for  the 
specific  time  periods  due  to  the  aperiodic  increase  and  decrease  of 
the  two  mean  currents  relative  to  each  other  as  they  develop  both  in 
depth  and  offshore  extent. 

Following  Renaud  (1986)  and  Batteen  et  al.  (1988),  the  damping 
(e-folding)  times  caused  by  model  heat  and  momentum  diffusion  was  also 
examined  to  see  how  they  compared  with  the  e-folding  times  from  the 
stability  analysis  of  the  mean  flow.  If  the  damping  time  was  of  the 
same  order  as  the  growth  rate  of  the  unstable  wave,  instability  of  the 
mean  flow  could  be  inhibited. 

Both  biharmonic  heat  and  momentum,  as  described  in  Section  II. A. 6, 
were  used  in  the  FE  model.  Following  Holland  and  Batteen  (1986),  the 
QG  baroclinic  mode  vorticity  and  thermal  equation  had  the  form: 

(V^  -  A^)^^  =  ...-A  V^V  +  A^B  vV.  (3-11) 
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Figure  3.29 


Baroclinic  instability  growth  rate  for  the 
unstable  jet  (from  Batteen  et  aJL.  ,  1988). 
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Figure  3.30 


Baroclinic  instability  growth  rate  for  time- 
averaged  periods  of  experiment  1. 
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where  V  was  QG  "temperature"  (rp^      -     V3  i"  two  layer  model),  A  was 
biharmonic  eddy  viscosity,   B  was  biharmonic  eddy  diffusion.  A  was 
inverse  Rossby  radius  and  only  the  damping  terms  were  kept  in  eqn. 
(3-10). 

Assuming  wave  numbers  (k,  Jl)  ,    (3-10)  became: 


V-t  =  -  [        1,2  .  p2  ,  .2         J  V".  (3-11) 


or       Vt  =  rp^e-r^.  (3-12) 


Therefore  the  e-folding  (damping)  time  of  baroclinic  modes  was: 
-1    (k^  +   J>^    +   X^) 

"'  -    A(k2  +     ^2)3  ^  ;^2g^i^2  ^  ^2^2  (3-13) 


The  damping  rates  (7"^)  and  time  scales  7  for  various  wavelengths 
using  values  of  A=B=2xlO^  ^  cm'*/sec,  as  in  the  PE  model  are  shown  in 
Table  5.  The  choice  of  values  A  and  B  are  consistent  with 
considerations  given  by  Semtner  and  Mintz  (1977),  who  found  that,  by 
using  a  biharmonic  operator,  an  optimum  value  to  use  can  be  selected 
which  will  both  control  computational  noise  and  reduce  the  diffusion 
of  mesoscale  features,  i.e.,  unrealistic  damping  will  be  reduced, 
except  for  features  which  are  poorly  resolved  while  the  advective 
effects  of  large-scale  features  will  be  enhanced.  They  used  the 
following  equation  to  obtain  damping  of  the  smallest  resolvable  wave 
in  either  coordinate  system  (e.g.,  when  k=7T/Ax  and  i=0): 

B  =  J  A  (Ax)^  (3-lA) 

where  B  was  the  biharmonic  diffusion,  A  was  the  Laplacian  diffusion 
and  Ax  was  the  grid  length.    Investigation  of  the  Reynolds  numl^er 
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yielded  the  results  that  A  cannot  be  less  than  1  x  10^  Gill's/sec  or  the 
Reynolds  number  will  be  too  large  and  unable  to  successfully  damp 
vorticity  on  a  grid  scale  of  20  km  (2Ax).  This  value  of  A  produces  e- 
folding  damping  times  that  are  too  short  (Table  5)  for  observing  the 
generation  of  reasonable  sized  mesoscale  features. 

From  Table  5,  the  damping  time  scale  of  ~22  days  correlated  with 
the  wavelength  (~70  km)  of  the  most  unstable  wave  in  the  imposed 
baroclinic  jet  case  of  Batteen  §t_  al.  (1988);  however,  the  baroclinic 
growth  rate  of  2.05  days  for  the  equivalent  wavelength  was  much 
shorter.  As  a  result,  the  diffusive  damping  should  not  have 
suppressed  any  development  of  instability.  Similar  results  were 
observed  for  the  longer  wavelengths  with  much  shorter  e -folding  time.s 
due  to  baroclinic  instability  than  diffusive  damping  times. 

Experiment  1  showed  different  results.  We  used  the  biharmonic 
heat  and  diffusion  to  damp  the  smaller  scale  "noise"  that  may  have 
masked  the  larger  scale  eddies  we  wanted  to  observe;  however,  in  this 
case  the  biharmonic  damping  may  have  inhibited  the  rapid,  baroclinic 
growth  of  more  intense  eddies.  The  most  unstable  wavelength  of  time- 
averaged  days  1-10  for  experiment  1  was  ~50  km  with  an  associated 
baroclinic  instability  growth  rate  of  ~4 . 7  days.  This  is  close  to  the 
biharmonic  damping  time  scale  of  5.8  days.  The  wavelength  of  ~44  km 
for  the  time-averaged  days  30-40  is  even  more  significant  with  a 
growth  rate  of  ~3.45  days  versus  the  damping  rate  of  ~3 . 5  days. 

Comparing  the  baroclinic  instability  growth  rates  (Figures  3.29 
and  3.30)  and  the  damping  rates  (Table  5),  we  can  conclude  some 
interesting  points  about  using  biharmonic  diffusion.   In  the  case  of 
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TABLE  5 
DAMPING  RATES  FOR  VARIOUS  WAVELENGTHS 


B I HARMONIC         LAPLACIAN 
(IxlO^^cmVsec)     ( IxlO^cm Vsec  ) 
L(km)         7  (days'M      j'^     (days)       7"^  (days)  7"'  (days) 


10  1.11x10^  .009  .02  .15 

20  0.67x10^  .15  .30  .58 

30  1.33x10^  .75  1.5  1.3 

40  4.17x10-^  2.4  4.7  2.4 


PE  MODEL 

PE  MODEL 

7  (days'M 

7"!  (days) 

1.11x10^ 

.009 

0.67x10^ 

.15 

1.33x10^ 

.75 

4.17x10"^ 

2.4 

1.72x10"! 

5.8 

8.33x10-2 

12 

4.55xl0"2 

22 

2.63x10-2 

38 

1.64x10-2 

61 

1.08x10-2 

93 

6.76x10-" 

1480 

50  1.72x10-'  5.8  12  3.7 

60  8.33x10-2  12  24  5.3 

70  4.55x10-2  22  45  7.2 

80  2.63x10-2  3g  7g  g  7 

90  1.64x10-2  61  122  12 

100  1.08x10-2  93  3^33  -^5 

200  6.76x10""  1480  2960  58 
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the  unstable  jet,  the  biharmonic  diffusion  damped  out  all  developing 
waves  ~50  km  or  less.  Experiment  1  also  had  all  developing  waves  50 
km  or  less  damped.  This,  unlike  Batteen  e_t  al.  (1988),  even  included 
the  most  unstable  wavelengths  (~40  km) .  With  the  shorter  wavelengths 
damped  out  or  growth  inhibited,  only  the  longer  waves  (~100  km)  with 
relatively  much  longer  damping  times  will  be  observed.  In  addition, 
the  time  when  they  actually  are  observed  are  even  longer  than  the  e- 
folding  baroclinic  growth  time  due  to  the  required  model  "spin-up" 
involved  in  the  development  of  the  mean  currents. 

Unlike  Laplacian  diffusion,  which  can  damp  all  wavelengths 
considerably  as  noted  by  Holland  (1978),  biharmonic  diffusion  should 
damp  the  shorter  wavelengths  (of  10-30  km)  in  a  short  amount  of  time 
and  the  longer  wavelengths  (longer  than  40  km)  in  a  longer  period  of 
time.  As  a  result,  ocean  models  with  very  fine  horizontal  resolution 
will  probably  have  shorter  wavelength  unstable  waves  damped.  Thus 
even  with  biharmonic  diffusion,  there  is  a  trade-off  between 
resolution,  damping  times  and  wavelengths,  and  for  very  fine 
resolution  (due  to  the  short  damping  times)  biharmonic  damping  may 
cause  shorter  wavelengths  which  are  resolvable  to  be  damped. 

Experiment  3,  which  included  the  wind  with  y-variation  in  wind 
stress  on  an  f -plane,  had  similar  stability  results  as  in  experiment 
1  due  to  the  interaction  of  the  coastal  jet  and  undercurrent. 
Experiment  2  with  wind  stress  curl  incorporated  did  not  go 
baroclinically  unstable  due  to  the  lack  of  a  developed  coastal  jet, 
which  as  showed  by  Batteen  et.  al^.  (1988),  was  necessary  for  baroclinic 
instability  to  occur.   A  strong  poleward  undercurrent  was  generated 
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that  eventually  reached  the  surface;  however,  the  vertical  shear 
associated  with  baroclinic  instability  was  not  present  so  that  no 
eddies  occurred.  The  results  of  experiment  2  paralleled  the  results 
of  Renaud  (1986),  who  generated  a  similar,  internal,  poleward  current 
that  likewise  did  not  develop  any  instability  (despite  satisfying  tlie 
necessary  conditions  for  both  barotropic  and  baroclinic  instability). 

The  last  experiments  (4,  5  and  6)  were  all  conducted  on  a  /3- plane. 
Holton  (1979),  based  on  a  zonal  flow  analysis,  investigated  the  effect 
of  /3  on  the  flow.  He  found  that  the  Beta  effect  strongly  stabilizes 
the  long  wavelengths  and,  as  expected,  the  flow  was  always  stable  for 
wavelengths  shorter  than  the  critical  wavelength,  L^,  .  The  long  wave 
stabilization  that  was  associated  with  the  Beta  effect,  was  a  result 
of  the  rapid  westward  propagation  of  packets  of  long  Rossby  waves. 

Olivier  (1987)  demonstrated  that  there  is  a  difference  in  flow 
behavior  between  a  non-zonal  and  zonal  flow.  He  found  that  for  a 
meridional  flow,  as  in  this  study  (representative  of  the  CCS),  energy 
can  be  released  without  any  component  of  /3  acting  on  it;  therefore  any 
shear  above  the  dissipation  level  may  produce  instability. 

Basic  instability  did  occur  in  the  /9-plane  experiments  4,  5  and  6, 
which  each  had  a  coastal  jet  overlying  a  poleward  undercurrent  with 
the  subsequent  development  of  eddies  and  jets.  These  and  the  previous 
experiments  provide  evidence  that  the  generation  of  complex  eddy  and 
jet  patterns  could  be  attributed  to  the  instability  created  by  the 
shear  between  the  coastal  jet  and  the  poleward  undercurrent. 
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IV.   COMPARISON  OF  MODEL  RESULTS  WITH  OBSERVATIONS 

A.   BRIEF  BACKGROUND  ON  THE  CCS 

Buyer  (1983)  described  the  classical  features  of  the  CCS  as 
consisting  of  a  baroclinic  alongshore  coastal  upwelling  jet  with  the 
strongest  equatorward  velocities  at  the  surface  over  the  midshelf  or 
outer  shelf  and  a  poleward  undercurrent  over  the  shelf  break.  Kosro 
(1986),  Kosro  (1987),  Huyer  and  Kosro  (1987)  and  Lynn  and  Simpson 
(1987),  however,  found  that  the  instantaneous  near-surface  currents 
often  deviated  substantially  from  the  time-averaged,  classical 
currents.  In  particular,  Kosro  (1987)  examined  synoptic  maps  of  the 
coastal  current  field  off  Northern  California  during  CODE  (Coastal 
Ocean  Dynamics  Experiment)  and  found  a  qualitative  correlation  between 
complex  temperature  patterns  in  satellite  imagery  and  intense  current 
structures  such  as  squirts,  eddies  and  jets.  Davis  (1985) 
investigated  CODE  drifter  buoy  results  and  concluded  that  it  was 
difficult  to  think  of  the  California  coastal  circulation  as  a  simple 
wind -driven  alongshore  current  with  cross-shelf  Ekman-driven 
circulation;  instead,  he  found  that  various  mesoscale  motions  could  be 
primary  features  for  cross -shelf  transport.  Other  observations  have 
shown  that  highly  energetic,  mesoscale  eddies  and  meandering  jets  can 
be  superimposed  on  the  broad,  slow  climatological  mean  flow  in  the  CCS 
(Bernstein  et  al..  ,  1977;  Mooers  and  Robinson,  1984;  Rienecker  et.  al .  . 
1985,  1988).    The  location  of  these  features  from  shore  and  their 
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intensity  can  also  vary  greatly,  depending  on  the  author  and  location 
of  the  observation  along  the  west  coast  of  the  United  States. 

B.   COMPARISONS 

A  comparison  of  model  results  with  available  observations  was 
carried  out  to  see  if  both  time -averaged  and  instantaneous  model 
simulations  of  the  coastal  jet,  undercurrent  and  eddies  were 
consistent  with  available  observational  data.  The  time -averaged  (over 
days  30  to  40)  comparisons,  prior  to  the  generation  of  eddies,  are 
shown  in  Table  6,  while  the  instantaneous  comparisons  to  highlight 
specific  characteristics  of  the  currents  and  eddies  are  shown  in  Table 
7.  The  observations  used  in  both  Tables  6  and  7  are  from  Huyer  and 
Kosro  (1987),  denoted  by  HK  in  column  1.  Table  6  also  includes  the 
seasonal  observations  from  Ikeda  and  Emery  (1984),  denoted  by  IE  in 
column  2.  The  time -averaged  results  of  Huyer  and  Kosro  (1987)  was 
obtained  from  a  set  of  synoptic  data  during  CODE  that  included  both 
strong  wind  events  and  relaxations.  These  observations  may  or  may  not 
be  representative  of  mean  climatological  conditions  in  the  CCS, 
therefore  the  comparison  between  the  observed  model  results  (using 
climatological  winds)  and  these  time-averaged  CCS  observations  can  be 
different.  In  this  study  the  wind  forcing  was  steady  with  either  an 
offshore  or  alongshore  variation  to  systematically  investigate  forcing 
of  the  model,  however  the  actual  ocean  regime  is  subject  to  both  the 
offshore  and  alongshore  variation  and  could  produce  CCS  currents  with 
slightly  different  characteristics.  The  letter  designations  used  in 
the  tables  have  the  following  meanings: 
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A.  Maximum  coastal  jet  velocity  (cms"^) 

B.  Offshore  location  of  coastal  jet  (km) 

C.  Offshore  extent  of  coastal  jet  (km) 

D.  Depth  of  inshore  coastal  jet  (m) 

E.  Maximum  vertical  shear  of  coastal  jet  (xlO"^  s"^) 

F.  Maximum  undercurrent  velocity  (cms'M 

G.  Offshore  location  of  undercurrent  axis  (km) 
H.   Width  of  undercurrent  (km) 

I .   Depth  of  undercurrent  axis  (m) 

J.   Maximum  zonal  eddy  diameter  (km) 

K.  Maximum  zonal  eddy  velocity  (cms"^). 
Most  of  the  time -averaged  model  results  of  experiment  1  (wind 
stress  without  curl  and  on  an  f -plane),  as  shown  in  Table  6.  compare 
quite  favorably  with  CCS  observations.  The  only  discrepancies  are  the 
following.  The  coastal  jet  is  slightly  deeper  and  weaker  than  the 
observations,  and  the  poleward  undercurrent  velocity  is  -10  cm/sec 
weaker  with  its  axis  location  ~10  km  closer  to  shore. 

The  utilization  of  a  flat-bottom  in  the  model  for  all  experiments 
versus  a  topography  including  a  shelf  and  slope  could  have  affected 
these  model  results.  McCreary  et.  aj^.  (1987)  found  that  the  presence 
of  a  shelf  tends  to  strengthen  the  coastal  jet  and  weaken  the 
undercurrent.  In  addition,  transient  rather  than  steady  wind  forcing 
could  result  in  a  more  realistic  undercurrent  (McCreary,  1987). 
Moreover,  our  value  for  the  average  alongshore  wind  stress  for 
experiment  1,  using  data  from  Nelson  (1977)  was  ~1  dyne/cm^,  which  was 
~33%  lower  than  the  calculated  values  used  by  Huyer  and  Kosro  (1987). 
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This  lower  value  for  wind  stress  would  also  contribute  to  a  weaker 
undercurrent  than  what  Huyer  and  Kosro  (1987)  observed. 

The  instantaneous  model  results  comparison  of  experiment  1,  as 
shown  in  Table  6,  also  shows  good  agreement  with  CCS  observations. 
Although  the  simulated  coastal  jet  is  still  deeper  and  weaker  compared 
with  observations,  the  model  poleward  undercurrent  is  more  corisistent 
with  the  instantaneous  observations.  The  largest  difference  between 
experiment  1  model  results  and  CCS  observations  is  in  the  maximum 
zonal  eddy  velocities.  CCS  observations  show  eddy  velocities  of  ~50 
cm/sec  while  the  model  shows  maximum  velocities  of  ~15  cm/sec.  This 
disagreement  could  likely  be  attributed  to  the  difference  between  the 
seasonal  wind  forcing  of  the  model  and  the  transient  event  wind 
forcing  observed  during  CODE. 

All  of  the  model  experiment  results  which  did  not  include  wind 
stress  with  curl  are  also  shown  in  Tables  6  and  7.  As  expected,  each 
of  the  experiment  results  compare  favorably  with  observations,  since 
overall  the  values  for  each  comparison  are  consistent  with  those  of 
experiment  1 . 

Since  the  wind  with  curl  experiments,  i.e.,  2  and  5,  did  not 
develop  prominent  surface  coastal  jets,  they  were  not  included  in 
Tables  6  and  7;  however,  experiment  5  did  show  evidence  of  a  very  weak 
(~2  cm/s)  coastal  jet  within  5  km  of  the  coast  near  the  end  of  the 
model  simulation  time.  The  dominant  current  features  of  both  the 
time-averaged  and  instantaneous  meridional  velocity  fields  were  a 
surface,  poleward  current  similar  to  the  Davidson  Current  and  a  broad 
surface,   equatorward  flow  ~100  km  offshore.    These   results  are 
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consistent  with  the  flat-bottomed,  wind  stress  with  curl  forcing  model 
results  of  McCreary  et  al.  (1987). 

Even  though  the  model  results  compare  favorably  both  with 
available  observations  (such  as  Huyer  and  Kosro,  1987)  and  with  other 
model  results  (McCreary  et.  al-  ,  1987),  it  is  difficult  to  establish 
which  experiment  is  the  best  simulation  of  the  CCS.  Based  on  these 
results,  probably  the  best  simulation  of  the  CCS  would  incorporate 
both  the  climatological  offshore  and  alongshore  variation  of  wind 
stress,  but  differences  may  still  occur  if  the  wind  stress  is  assumed 
to  be  time  -  independent .  Observations  show  that  the  coastal  jet, 
undercurrent,  eddies  and  Davidson  Current  can  vary  greatly  both 
spatially  and  temporally.  A  model  that  can  correctly  simulate  the  CCS 
and  its  complex  features  would  require  accurate  observational  data, 
parameterizations ,  boundary  conditions  and  a  host  of  other  conditions 
too  numerous  and  complex  for  the  scope  of  this  study. 
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V.   SUMMARY  AND  RECOMMENDATIONS 

A .   SUMMARY 

This  study  used  a  high-resolution,  multi-level  PE  ocean  model  to 
investigate  wind  forcing  as  a  possible  generation  mechanism  for 
mesoscale  eddies  and  jets  in  the  CCS.  A  band  of  seasonal,  steady 
winds,  either  with  or  without  wind  stress  curl,  and  either  with  or 
without  alongshore  variability,  was  used  as  model  forcing  on  either  an 
f -plane  or  a  ;0-plane  in  an  idealized,  flat-bottomed  oceanic  regime 
along  an  eastern  ocean  boundary.  In  addition,  a  stability  analysis 
was  made  to  determine  if  both  the  necessary  and  sufficient  conditions 
for  instability  processes  to  occur  were  satisfied.  The  analysis 
showed  that  the  necessary  condition  for  barotropic  instability  and 
both  the  necessary  and  sufficient  conditions  for  baroclinic 
instability  were  satisfied  in  some  of  the  experiments. 

The  model  results  of  experiment  1,  which  included  wind  stress 
without  curl  on  an  f -plane,  showed  the  development  of  an  equatorward 
coastal  jet  and  poleward  undercurrent.  Baroclinic  instability  due  to 
the  unstable  jet  and  undercurrent  occurred  after  ~^0  to  45  days 
resulting  in  the  generation  of  cyclonic  and  anticyclonic  eddies  and 
jets.  Similar  results  occurred  in  experiment  4,  which  had  the  same 
form  of  wind  stress  as  in  experiment  1,  but  used  a  /9-plane  rather  than 
an  f -plane.  Due  to  the  Beta  effect,  Rossby  waves  propagated  offshore 
as  a  packet,  with  longer  followed  by  shorter  waves  moving  westward. 
Due  to  the  use  of  vertical  mixing  in  the  model,   consistent  with 
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McCreary  e^  al.  (1987),  the  Rossby  waves  did  not  advect  either  the 
coastal  jet  or  eddies  offshore. 

Experiment  2,  which  included  wind  stress  with  curl  on  an  f -plane, 
showed  initially  the  development  of  a  coastal  jet.  The  implementation 
of  a  steady,  positive  wind  stress  curl  near  shore  resulted  in  the 
generation  of  both  a  deep,  broad  poleward  surface  current  near  the 
coast,  similar  to  the  Davidson  Current,  and  an  equatorward  surface 
flow  located  farther  offshore.  No  eddies  or  jets  developed  in 
experiment  2  due  to  the  lack  of  an  unstable  vertical  shear  between  a 
coastal  surface  jet  (which  was  replaced  by  the  poleward  undercurrent) 
overlying  a  poleward  undercurrent.  Experiment  5  also  incorporated 
wind  stress  curl,  but  used  a  ^-plane  rather  than  an  f -plane.  Due  to 
the  presence  of  the  ^-plane,  a  coastal  jet  inshore  of  the  Davidson 
Current  developed  late  in  the  model  simulation  time  and  within  this 
jet,  eddies  were  generated.  As  expected  with  the  Beta  effect,  a 
packet  of  Rossby  waves  propagated  offshore,  as  in  experiment  4. 

The  model  results  of  experiment  3.  which  had  a  y-variation  in  the 
alongshore  wind  stress  without  curl  on  an  f- plane,  were  comparable  to 
experiment  1.  The  variation  in  alongshore  wind  stress  plaved  a  key 
role  in  determining  the  location  of  the  eddy  development  region.  In 
particular,  eddies  were  generated  in  the  localized  area  of  maximum 
alongshore  wind  stress  just  south  of  the  Mendocino  region.  Experiment 
6  also  incorporated  a  y-variation  in  the  alongshore  wind  stress,  but 
used  a  /3-plane  rather  than  an  f -plane.  The  fact  that  eddies  developed 
further  north  of  the  localized  eddy  generation  area  of  experiment  3 
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showed  that  the  Beta  effect  can  also  play  a  key  role  in  determining 
the  location  of  eddy  generation  regions. 

A  comparison  of  model  results  with  available  observations  showed 
that  both  time -averaged  and  instantaneous  model  simulations  of  the 
coastal  jet,  undercurrent  and  eddies  were  consistent  with  available 
CCS  observations  (e.g.,  Huyer  and  Kosro,  1987)  and  other  model  results 
(McCreary,  et.  ai- .  1987).  The  instantaneous  model  simulations  of  the 
CCS  duplicated  the  great  variability  in  time  and  space  of  the 
observed  CCS  while  the  time-averaged,  simulations  showed  a  classical 
two-dimensional  coastal  jet. 

B.   RECOMMENDATIONS 

The  results  from  these  experiments  strongly  support  the  hypothesis 
that  wind  forcing  can  be  a  significant  generation  mechanism  for  eddies 
and  jets.  It  should  be  noted,  however,  that  this  study  employed  the 
constraints  of  a  regular,  straight  coastline  and  a  flat-bottom.  It 
also  made  use  of  steady,  idealized  wind  stress  patterns.  Batteen  e_t 
al .  (1988),  using  the  same  PE  model,  imposed  a  climatological  mean 
coastal  jet  and  undercurrent  over  an  idealized  Mendocino  Escarpment 
and  Ridge  complex,  and  showed  that  topography  can  influence  the 
characteristics  of  eddies.  McCreary  et.  al-  (1987),  found  that  the 
presence  of  a  coastal  shelf  in  an  ocean  model  can  strengthen  the 
coastal  jet  and  weaken  the  undercurrent.  Ikeda  and  Emery  (1984) 
concluded  that  current  meanders  could  be  triggered  by  alongshore 
variations  (capes)  in  the  coastline  of  California  and  Oregon,  and  grow 
due  to  baroclinic  instability  of  the  mean  currents.    Based  on  the 
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results  from  these  studies,  future  studies  should  include  both  an 
irregular  coastline  and  bottom  topography. 

In  addition,  different  open  boundary  conditions  as  discussed  by 
R<^ed  and  Cooper  (1985),  should  be  incorporated  and  tested  in  the 
model.  In  particular,  wind  forcing  with  both  a  local  solution  and  a 
global  solution  (R(^ed  and  Smedstad,  1984)  could  be  implemented  and 
tested.  This  latter  incorporation  may  permit  the  wind  forcing  to  be 
applied  over  the  entire  domain  (even  at  the  open  boundaries),  while 
still  allowing  the  free  propagation  of  waves  from  the  boundaries 
througli  the  computational  domain. 

Time-dependent  wind  forcing,  such  as  wind  events  and  relaxations, 
should  also  be  systematically  investigated  to  see  if  transient  wind 
forcing  can  be  an  important  eddy  generation  mechanism.  McCreary  et 
al ■  (1987)  has  shown  that  the  inclusion  of  both  annually  periodic  and 
remote  winds  can  significantly  influence  the  mean  flow  currents  of  the 
CCS.  Future  studies  therefore  should  include  time -dependent  winds 
along  with  remote  forcing. 

Future  experiments  should  also  include  not  only  the  alongshore 
variation  in  wind  stress,  but  also  the  offshore  variation  in  wind 
stress.  The  recent  higher  resolution  wind  data  of  Bakun  (1988)  could 
be  incorporated  as  the  wind  forcing  and  systematically  investigated. 
This  should  result  in  even  more  favorable  model  comparisons  with 
available  CCS  observations. 
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